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T-H Abstract 

One place where Belyi functions occur is puUback transformations 
of hypergeometric differential equations to Fuchsian equations with few 
^ singularities. This paper presents a complete classification of Belyi 

Q functions for transforming certain hypergeometric equations to Heun 
equations (i.e., canonical Fuchsian equations with 4 singularities). The 
\^ considered hypergeometric equations have the local exponent differences 

T-H l/k,l/£,l/m that satisfy k,£,m G N and the hyperbolic condition 

1/k + l/£ + 1/m < 1. In total, we find 872 such Belyi functions up 

O' ' to Mobius transformations, in 366 Galois orbits. Their maximal degree is 

60, which is well beyond reach of standard computational methods. To 
obtain these Belyi functions, we developed two efficient algorithms that 
T exploit the implied hypergeometric-to-Heun transformations. 

1 Introduction 



> 



Belyi functions and dessins d'enfants [241 is a captivating field of research 
in algebraic geometry, complex analysis, Galois theory and related fields. 
However, computation of Belyi functions of degree over 20 is still considered 
00 hard [TH Example 2.4.10] even for genus Belyi coverings — > P^. This 

CO computational difficulty promises a long lasting appeal, both for constructivists 

psj and theoreticians. Grothendieck [HI pg. 248] doubted that "there is a uniform 

T-H method for solving the problem by computer". The subject of this paper is 

effective computation of certain Belyi functions V^, of degree up to 60, 

I utilizing the fact that those functions transform Fuchsian differential equations 

^ with a small number of singularities. 

k> In this paper, a Belyi functioi^ (or a Belyi covering) is a rational function 

if : P], ^ Fl that branches only in the 3 fibers z = ip{x) e {0,1, oo}. We 
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^Generally, a Belyi function is a map tp : S ¥]. from a Riemann surface S to P^(C) that 
only branches in the fibers z £ {0, 1, oo}. Definitions 1 1 . l[|l . 3| could be applied to general Belyi 
maps S — > P^. Often a Belyi function is defined by requiring branching in at most 3 fibers 
rather than in the fixed fiber set {0, 1, oo}. Although any 3 fibers can be m oved to {0, 1, oo} 



by a Mobius transformation of P , the definition difference is significant in {4.5 
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distinguish the two projective Hnes by the indices x, z just as in [5^. To describe 
the Belyi functions we classify, we introduce the fohowing definitions. 

Definition 1.1. Given positive integers k,£,m, a Belyi function (p : V]. ^ V], 
is called {k, £, m)-regular if all points above z — 1 have branching order k, all 
points above z — have branching order £, and all points above z = oo have 
branching order m. 

Examples of (2, 3, m)-regular Belyi functions with m G {3, 4, 5} are the well- 
known Galois coverings of degree 12, 24, 60 with the tetrahedral 
A4, octahedral ^4 or icosahedral ^5 monodromy groups, respectively. Platonic 
solids give their dessins d'enfant [TS] . 

Definition 1.2. Given yet another positive integer n, a Belyi function : P!^ — >■ 
V], is called (/c, ^, m)-minus-rt- regular if with exactly n exceptions, all points 
above z = 1 have branching order k, all points above z = have branching order 
£, and all points above z = 00 have branching order m. We will also use the 
shorter term (fc, £, m)-minus-n. 

Examples of (k, i, r7i)-minus-2 functions are quotients of the just mentioned 
Galois coverings by a cyclic monodromy group. If l/k + l/£+l/m > 1 and n > 3, 
there are (fc, £, m)-minus-n Belyi functions of arbitrary high degree. They give 
Kleinian pull-back transformations [TOl HH] to second order Fuchsian equations 
with finite monodromy (i.e., a basis of algebraic solutions) from a few standard 
hypergeometric equations. An example of a (2, 3, 5)-minus-3 Belyi function of 



degree 1001 is given online at [T71 NamingConvention]. As Remark 5.1 here 



shows, (fc, m)-minus-l-regular Belyi functions exist only if 1 G {k,£,m}. 

Definition 1.3. A Belyi function Lp is called minus-n-hyperbolic if: 

(i) there are positive integers k,£,m satisfying 1/k + l/£ + 1/m < 1 (the 
hyperbolic condition^ such that ip is {k, £,m)-minus-n-regular; 

(ii) there is at least one point of branching order k above z = I, a point of 
order £ above z = 0, and a point of order m above z = 00. 

Minus-3-hyperbolic Belyi functions are listed in [20l §9] . Table 3 in [20^ lists 
nincj^ Galois orbits of such Belyi functions, of degree 6, 8-10, 12, 18 or 24. 

This paper gives all minus-A-hyperbolic Belyi functions P^ — ^ P^. The 
motivation is that they give transformations of Gauss hypergeometric dif- 
ferential equations without Liouvillian ^24J solutions to Heun equations (i.e., 
Fuchsian equations with 4 singularities). This allows to express non-Liouvillian 
Heun functions in terms of better understood Gauss hypergeometric functions. 
The application to these transformations of Fuchsian equations is discussed in 
^ This paper, combined with the list of parametric hypergeometric-to-Heun 
transformations in |23j . covers all non-Liouvillian cases of hypergeometric-to- 
Heun transformations. 

We used two algorithms to compute the minus-4-hyperbolic Belyi functions. 
They both utilize the fact that these Belyi functions give hypergeometric-to- 
Heun transformations. One algorithm is probabilistic and uses modular lifting. 



■^Minus-3-hyperbolic Belyi functions give rise to the hypergeometric transformations 
described in |20l §9], There are 10 different such Belyi functions up to Mobius transformations, 
in 9 Galois orbits. The degree 18 Belyi function there is defined over 7). 



2 



It exploits the fact that Heun's equation is represented by few parameters. The 
other algorithm is deterministic, and uses existence of a hypergeometric-to-Heun 
transformation to get more algebraic equations for the (a priori) undetermined 
coefficients of a Belyi function. 

The branching patterns are enumerated in ^|3j following the approach from 
[22]. Section [4] discusses special obstructed cases of encountered Belyi functions. 
The application to hypergeometric-to-Heun transformations is explained in ^|5j 
Our algorithms are presented in ^ The Appendix sections give ordered lists 
A-J of computed Belyi functions, discusses composite Belyi functions, and 
compares our results with Felixon's list of Coxeter decompositions in the 
hyperbolic plane. All dessins d'enfant of computed Belyi coverings are depicted 
in this paper, most of them next to the A~J tables of ^B] Our list of dessins is 
larger than [TJ HJ [T3] combined. 



2 Organizing definitions, examples 

We start with a few definitions that will help us to organize the list of Belyi 
functions. Then we take a relaxed look at a few examples, including those of 
the largest degree 60. At the same time, dessins d'enfant are introduced more 
fully, setting a geometric tone of our presentation. 

Definition 2.1. Let ip be a {k,£^m)-minus-n-regular Belyi function for some 
n. The regular branchings of ip are the points above z = 1 of order k, the points 
above z ^ of order i, and the points above z — 00 of order m. The other n 
points in the three fibers are called exceptional points of tp. A branching fraction 
of ip is a rational number A/B, where A is a branching order at an exceptional 
point Q, and B G {k,i,m} is the prescribed branching order for the fiber of Q. 

Definition 2.2. Let ip : Vl ^ be a {k,£,m) -minus-A Belyi function. Let 
9i J <Z2 , 93 7 94 denote its exceptional points. T/ie j-invariant of is the j -invariant 



of the elliptic curve — Yiqi^ooi-^ ~ g'i'vcn by formula (2.2) below. It is 
invariant under Mobius transformations of V], . 

A canonical form of ip is a composition of ip with a Mobius transformation 
that has three exceptional points at x = 0,1,00. The fourth exceptional point 
then becomes x = t, where t is a cross-ratio of qi,q2,q3,q4. The cross-ratio 
depends on the order 0/91,92,93,947 Cind there is an S^-orbit (S3 = S4/V4) 

< 1 1 11 , , 

t,l-t, ,-, -1--^ 2.1 

t-l i' 1-t tj ^ ' 

of related cross-ratios. Any of these values is a i-value of ip. The j-invariant is 

256(f^-t + l)3 

As an example, t e {-1, 2, 5} gives j ^ 1728. If j ^ {0, 1728} then the six 
values in above S'3-orbit are distinct. The S'3-action gives a homomorphism 
S3 ^ Gal(Q(0/Q(j))7 hence [Q(i) : q{j)] e {172,6}. 

Definition 2.3. The t-field resp. j-ficld of ip is the number field generated by 
a t-value resp. the j-invariant. The r -field (canonical realization fieldj of ip is 
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the smallest field over which a canonical form of if is defined. These fields do 
not depend on the ordering of the 4 exceptional points. 



Example 2.4. The degree 12 rational function 

Mx^{x-if{x-<d) 



27{x - l)(8a;3 - 72x^ ~ 27x + 27f 



is a (2, 3, 9)-mmMS-4 Belyi map. Indeed, the numerator of ip{x) — 1 is a full 
square. It is already in a canonical form, as x = 0, x = 1 and x = od (of 
branching order 2) are exceptional points. The fourth exceptional point x ~ 9 is 



a t-value. The j-invariant is equal to 2^73^/3* by formula (2.2). 

The branching pattern of ip is given by three partitions of the degree d — 12 
into branching orders above 0,1, cx). Using the notation in [22\, we express the 
branching pattern of ip shortly as follows: 

6 [2] = 3 [3] + 2 + 1 = [9] + 2 + 1. 

The prescribed branching orders are indicated with square brackets, with their 
multiplicity in front. The 4 branching orders that are not enclosed in square 
brackets represent the 4 exceptional points. Dividing them by their prescribed 
branching order(s) produces the 4 branching fractions: 1/3,2/3,1/9,2/9. 

In the appHcation setting of hypergeometric-to-Hcun transformations in ^|5] 
the regular branchings will become regular points (after a proper projective 
normalization) of the pulled-back Heun equation H; the exceptional points will 
be the singularities of H; and the branching fractions will be the exponent 
differences of H. The exponent differences of the hypergeometric equation under 
transformation will be 1/fc, 1/m. Example 2.4 will be continued in Sjsj 



Definitions |2.2[ |2.3| will be used to group the obtained Belyi functions into 
manageable classes. The Belyi functions will be listed twice in this paper. The 
first list is Tables 2.3.7-3.4.4 of Sjs] Its ordering by the (fc, m)-triples and 
branching patterns refiects the classification scheme. In Appendix ^B] the list of 
Galois orbits is grouped and ordered by the j-fields, t-fields, branching fractions. 
This order allows quick recognition whether a given Heun function is reducible 
to a hypergeometric function with a rational argument Lp. 

Belyi functions nicely correspond to certain graphs called dessins d'enfan^ 
Mimicking [4t Section 2] , we spell out standard correspondences for genus Belyi 
functions as follows. There are 1-1 correspondences between these objects: 

(I) Belyi functions V\ up to Mobius transformations x {ax + h) /{cx -f d). 



•^Generally, a dessin d'enfant 1241 is a bi-colored graph (possibly with multiple edges), with 
a cyclic order of edges around each vertex given. This defines a unique (up to homotopy) 
embedding of the bi-colored graph into a Riemann surface. Customarily, the vertex colors are 
black and white. The dessins d'enfant of genus Belyi coverings can be drawn on a plane, 
as the Riemann sphere minus a point is homeomorphic to a plane. Given a Belyi covering ip, 
its dessin d'enfant is realized as the pre-image of the interval segment [0, 1] C M C C onto its 
Rieman surface, with the vertices above z = colored black and the vertices 2 = 1 colored 
white. The branching pattern of ip determines the degrees (i.e., valencies) of vertices of both 
colors of its dessin d'enfant, and the degrees of cells on the Riemann surface. The cell degree 
is determined by counting vertices of one color while tracing its boundary. The degree of a 
dessin d'enfant is the degree of the corresponding Belyi function. 
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(II) Plane dessins d'enfant, up to homotopy on the Riemann sphere. 

(III) The triples {go, gi, goo) of elements in a symmetric group Sd, such that: 

• .9o.9iffoo = 1; 

• the total number of cycles in go,gi,goo is equal to d + 2; 

• ffOiffiiSoo generate a transitive action on a set of d elements; 

up to simultaneous conjugacy of go,gi,goo hi Sd- 

(IV) Field extensions of Q_{z) of genus 0, unramified outside z — 0,1, oo. 

Part (III) gives the monodromy presentation of a Belyi covering, and d is 
the degree. The dessins d'enfant is basically a graphical representation of the 
combinatorial data in (III) . This paper presents all obtained dessins pictorially, 
while the accompanying website [17] gives the Belyi maps (I) , the permutations 
in (III) and other data (such as j, r-fields). The number c? + 2 is illuminated 



in the proof of Lemma 3.1 For each fiber z g {0, 1, oo}, the conjugacy class of 
gz is determined by the partition of d that reflects the branching pattern in the 
fiber. Part (IV) is convenient for considering the composition structure of Belyi 
maps; see Appendix [Cj 

The considered Belyi functions have rather regular dessins d'enfant. Deflni- 
tions |l.l]|1.3| are easy to reformulate for dessins d'enfant; 

Definition 2.5. A dessin d'enfant is called (fc, TO)-minus-n-regular if, with 
exactly n exceptions, all white vertices have degree k, all black vertices have 
degree £, and all cells have degree m. 

Definition 2.6. A dessin d'enfant T is ca/Zed minus-n-hyperbolic if: 

(i) there are positive integers k,£,m satisfying l/k + l/£ + l/m < 1 such that 
r is (k,£,m)-minus-n-regular; 

(ii) there is at least one white vertex of degree k, a black vertex of degree £, 
and a cell of degree m. 

All minus-4-hyperbolic dessins d'enfants could be found by a combinatorial 
search on a computer. But with our Maple implementations it was faster to 
compute first the minus-4-hyperbolic Belyi functions, and then compute their 
monodromy permutations in (III). This paper presents all minus-4-hyperbolic 
dessins (up to complex conjugation), most of them next to the tables of 
Appendix ^ 

In total, there are 872 Belyi functions of the minus-4-hyperbolic type, up 
to Mobius transformations in both x and z. They come in 366 Galois orbitf]^ 
In leap years we could decorate a calendar with the minus-4-hyperbolic dessins 
d'enfant, one Galois orbit per day. We categorize and label the Galois orbits of 
the objects in (I)"(IV) as A1-J28; see j |3.1| and Appendix S|b] The largest Galois 
orbit J28 has 15 dessins, for a (2,3, 7)-minus-4 branching pattern of degree 37. 
Completeness is checked with two independent algorithms and other checks, see 
^and Appendix 

*Belyi functions are explicitly defined over algebraic number fields, and the absolute 
Galois group Gal(Q/Q) permutes Belyi coverings with the same branching pattern. The 
size of a Galois orbit of dessins d'enfant is the degree of the moduli field; see ^ Given 
a branching pattern, the set of Belyi coverings with that branching pattern is finite (up to 
Mobius transformations), possibly empty. The Galois action does not need to be transitive 
on this set, and several Galois orbits with the same branching pattern may appear. 
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Figure 1: The degree 60 dessins d'enfant 



The highest degree of a niinus-4-hyperbohc Belyi function is 60. Its 
branching pattern is 30 [2] = 20 [3] = 8 [7] + 1 + 1 + 1 + 1. There are two 
Galois orbits for this branching pattern, with three dessins each. We identify 
the two Galois orbits as H14 and H46. The dessins d'enfant for these Belyi 
functions are depictecj^ in Figure [ij The 4 exceptional points in each dessin are 
represented by circular loops; they could be assumed to lie in the center of each 
cell of degree 1. The other cells (including the outer ones) have degree 7. The 
left-most dessins of H14 and H46 clearly have a reflection symmetry, hence they 
are defined over M. The other two dessins of H46 are mirror images of each 
other, and are related by the complex conjugation. 

The Belyi functions of degree 60 are composite. H46 is a composition of a 
degree 15 Belyi function H45 with two consecutive quadratic substitutions. The 
intermediate Belyi function J19 of degree 30 has a Galois orbit of size 9, while 
H45 is defined over the same cubic field as H46. The branching fractions of H45 
are 1/2,1/2,1/2,1/7. The dessins d'enfant for the component functions H45 
and J19 are depicted in Figure [2] 

The covering H14 is a composition of degree 30 Belyi function HIO with a 
quadratic substitution. The functions HIO and J19 have the same branching 
fractions 1/2,1/2,1/7,1/7. The Belyi functions HIO, H14 are examples that 

^The dessins in Figure [l] have all white vertices of order 2, hence they are examples of 
clean dessins d'enfant. It is customary to depict clean dessins without white vertices, so that 
edges connect black vertices directly, and loops are possible. A white vertex is then implied 
in the middle of each edge. Our policy of drawing dessins is the following. White vertices 
of order 2 are not shown, but the edges going through them are drawn thick. Other white 
vertices are shown, but the incident edges are drawn thin. A black vertex of degree > 2 is 
not drawn (as it is a clear branching point), unless it is incident to a thin edge. Take a look 
at Figure [2] to see effects of this policy. Figures ^ |2] depict all dessins d'enfant for H14, H45, 
H46, J19, including complex-conjugated pairs. In other pictures, we will not draw dessins 
whose complex conjugates are already present. 
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H45 



J19 




Figure 2: The degree 15 and 30 components of H46. 

have an obstruction, as described in J|4] This has interesting geometric 
consequences for the dessins d'enfant. Although both have a totahy real moduli 
field Q(cos^), not all dessins of HIO and H14 have a reflection symmetry. 
Rather, the complex conjugation may give a homeomorphic dessin, identifiable 
with the original only after an automorphism of the Riemann sphere. For 
example, consider the middle and the right-most dessins of H14 in Figure [T] 
The dessins d'enfant for HIO are depicted in Figure [3] together with most of 
other examples with an obstruction. 



3 The branching patterns 

We enumerate the possible branching patterns in the same way as was done for 
parametric hypergeometric-to-Hcun transformations in 22]. To end up with a 
finite number of cases, we use Hurwitz formula and the hyperbolic condition 
l/k + l/£+ 1/m < 1. Without loss of generality, we assume the non-decreasing 
order k < £ < m for the regular branching orders from now on. 

Lemma 3.1. Let ip be a minus-A-hyperbolic Belyi covering of degree d, with the 
regular branching orders k < £ < m d Z>o- Then 

(i) There are exactly d ~ 2 regular branchings and 4 exceptional points. 
, . , d d d 

(''^'-[-k\~b\-[^\^'- 

2-S 



(Hi) Let S denote the sum of 4 branching fractions. Then d - 



1 
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(iv) (^^-^ - - + 4 < 0. 

/ > 1 1 1 , 

Proof. By Hurwitz formula (or PHI Lemma 2.5]), there are 3d — {2d — 2) = d + 2 
distinct points above {0, 1, oo} when is a Belyi map — P-'^. The first claim 
follows. The number of regular branchings is at most [d/fcj + ld/£\ + [d/mj. 
This gives the inequality in (ii). The number d — 2 of regular branchings is also 
equal to d/k + d/£ + d/m — S, giving the degree formula in (Hi). 



We have d > m, otherwise condition (ii) of Definition 1.2 is not satisfied. 
Combining this with the degree formula gives the inequality in (iv). Together 
with m > 4, the inequality in (iv) gives 1 — 1/€ < 1/2. Part (v) follows. □ 

The inequalities in (iv), (v) gives a finite list of triples (k, £, m). Setting S = 
4/m in part (Hi) gives an upper bound for d, leaving the following candidates 
for {k, t, TO, d): 

(2, 3, 7, < 60), (2, 3, 8, < 36), (2, 3, 9, < 28), (2, 3, 10, < 24), 
(2,3,11, < 21), (2,3,12,<20), (2, 3, 13, < 18), (2, 3, 14, < 18), 
(2, 3, 15, < 17), (2, 3, 16, < 16), (2, 4, 5, < 24), (2, 4, 6, < 16), 
(2,4,7,<13), (2,4,8,<12), (2,4,9,<11), (2, 4, 10, < 10), 
(2, 5, 5, < 12), (2, 5, 6, < 10), (2,5,7,<9), (2,5,8,<8), 
(2,6,6,<8), (2,6,7,<7), (3,3,4,<12), (3,3,5,<9), 
(3,3,6,<8), (3,3,7,<7), (3,4,4,<6), (3,4,5,<5), 
(4,4,4, < 4). 

The last two candidates give less than 4 exceptional points. 

Given a candidate tuple (fc,£, TO,d), it is straightforward to enumerate the 
correspondent branching patterns. Let hi,ho,hoo denote the eventual number 
of regular branchings in the fibers z — 0,l,cx), respectively. Then we have 
hi + Hq + hrxi =d— 2 by Lemma |3.1| (^zJ, and < hi < [d/k\, etc. With a 
possible integer solution {hi, /iq, hoc) at hand, we have to partition the numbers 
d— khi, d — Iho, d— mhao into total 4 positive parts, not equal to the respective 
regular orders fc, m. For example, if {k, £, m, d) = (2, 3, 7, 28) then we are lead 
to partition either 7 into 3 parts (in the to = 7 fiber) or 4 into 4 parts (in 
the £ ^ 3 fiber). There are four such partitions of 7 and one of 4, hence five 
branching patterns. 

In total, there are 378 branching patternsj^ for minus-4-hyperbolic Belyi 
functions. We list them in the first two columns of Tables 2.3.7-3.4.4, by giving 
their branching fractions and the degree. The table titles refer to the tuple 
{k,i,m). The branching fractions are left unsimplified (e.g. 4/8 instead of 1/2) 
to keep the fibers and branching orders of exceptional points plainly visible. 
The branching patterns are uniquely determined by the unsimplified branching 
fractions. 

^One branching pattern 6 [3] = 9 [2] = 8 + 7 + 1 + 1 + 1 is counted twice. It appears in 
Tables 2.3.7 and 2.3.8 because it is (2, 3, m)-minus-4-regular for two values of m (=7 or 8). It 
turns out, however, there are no Belyi functions with this branching pattern. 
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3.1 Summary of computed results 

The last two columns of Tables 2.3.7-3.4.4 give information about found 
coverings with the considered branching patterns. Computation of Belyi 
coverings for each possible branching pattern is the most demanding step. The 
algorithms used to generate and verify the list of Belyi functions are presented 
in Section fj6] 

The third column of Tables 2.3.7-3.4.4 gives a label for every Galois orbit 
with the branching pattern defined by a sequence of 4 branching fractions in 
the first column. With the application to Heun equations in mind, we group 
the Belyi functions by the Q-extension of the j-invariant. The cases with j G Q 
are further grouped by the i-field. We group the computed minus-4-hyperbolic 
Belyi functions into the following 10 classes, labeled alphabetically from A to J: 



A1-A24 
B1-B34 
C1-C42 
D1-D50 
E1-E25 

F1-F25 
G1-G52 
H1-H53 
11-133 
J1-J28 



the Belyi functions with j = 1728, that is t e {-1, 2, 1/2}; 
the other Belyi functions with t G Q; 

the Belyi functions with j E Q and a real quadratic i-field; 

the Belyi functions with j G Q and an imaginary quadratic t-field; 

the Belyi functions with j € Q, when the t-field is the splitting 

field (of degree 6 always) of a cubic polynomial in Q[z]; 

the Belyi functions with a real quadratic j-field; 

the Belyi functions with an imaginary quadratic j-field; 

the Belyi functions with a cubic j-field; 

the Belyi functions with a j-field of degree 4 or 5; 

the Belyi functions with a j-field of degree at least 6 (and < 15). 



In each class, the Galois orbits of Belyi functions are ordered by the criteria 
described in Appendix SjXj A numbered label refers to a whole Galois orbit of 
Belyi functions (or dessins d'enfant), as already mentioned in Sj2] If there is 
more than one Galois orbit with the same branching pattern, a line is devoted 
to each Galois orbit in Tables 2.3.7-3.4.4. 

The last column of the same tables gives basic information about the size 
of Galois orbits, j-fields, t-fields of the computed Belyi functions. The j-ficld is 
indicated as follows: 

• by the field degree d, in the power notation j'^; 

the fiek 
+ aX^ + bX + c- 



• if the field is quartic, j'^{a,b,c) conveys a minimal field polynomial 



• if the field is cubic, j (a, 6) conveys a minimal field polynomial X +aX+b; 

• if the field is quadratic, j^(-\/a) means the field Q(-y/a); 

• if j = 0, it is stated so; 

• for j e Q \ {0}, no j-notation is given, but the i-field and (possibly) the 
moduli field are indicated. 

The i-field is specified as follows: 

• if the j-field is indicated, the i-field degree d is given (in the power notation 
f^) only if j 7^ and the t-field is an extension of the j-field; 
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Table 2.3.7. 


1/7,1/7,1/7,1/7 60 


H14 


J^(-7,7),i6 




H46 


j3(-7,14),ti8 


^/l}/7.:^/7?/7 54 


D28 


f(V-3) 


1/3,1/7,1/7,1/7 52 


D35 




1/7,1/7,1/7,3/7 48 
1/7,1/7,2/7,2/7 


— 


no covering 


A21 


t = -1 




15 


/(-7,0,14),t« 


1/3,1/7,1/7,2/7 46 


J7 




1/2,1/7,1/7,1/7 45 
1/3,1/3,1/7,1/7 44 


G44 


J^(a/=7), 


Hll 


J^(-7,7),i6 




J26 


jl3, ^26 


1/7,1/7,1/7,4/7 42 


G33 


j'(y=7), 


1/7,1/7,2/7,3/7 


— 


no covering 


1/7,2/7,2/7,2/7 


— 


no covering 


1/3,1/7,1/7,3/7 40 


J16 




1/3,1/7,2/7,2/7 


G34 




1/2,1/7,1/7,2/7 39 


G43 


j2(y=7), f4 


1/3,1/3,1/7,2/7 38 


J13 




2/3,1/7.1/7.1/7 


16 


/(-7,0, 14), 


1/2,1/3,1/7,1/7 37 


J28 




1/3,1/3,1/3,1/7 36 


114 


/(0, 14,21),t24 


1/7,1/7,1/7,5/7 
1/7,1/7,2/7,4/7 


— 


no covering 


A22 


m2(v^),t=-l 




D40 


i(V=7) 


1/7,1/7,3/7,3/7 


B12 


t = 32 




128 




1/7,2/7,2/7,3/7 

2/7,2/7,2/7,2/7 


— 


no covering 


A23 


t = -1 




E22 


t^P'(9, 2) 


1/3,1/7,1/7,4/7 34 

1/3,1/7,2/7,3/7 

1/3,2/7,2/7,2/7 


D17 




JIG 






no covering 


1/2,1/7,1/7,3/7 33 


132 


^5 flO 
J , 


1/2,1/7,2/7,2/7 
1/3,1/3,1/7,3/7 32 


G41 


j2(y=7), ^4 


AlO 


i = -1 




F21 


.7'(^/7), ^4 


1/3,1/3,2/7,2/7 


C6 


t(V3) 

J , 




123 


2/3,1/7,1/7,2/7 


G32 




1/2,1/3,1/7,2/7 31 
1/2,1/2,1/7,1/7 30 


J24 




HIO 






J19 


f, ^18 


1/3,1/3,1/3,2/7 


— 


no covering 


1/3,2/3,1/7,1/7 


J6 




1/7,1/7,1/7,6/7 


D22 


J=0 


1/7,1/7,2/7,5/7 
1/7,1/7,3/7,4/7 


D5 


^(^/=l) 




no covering 


1/7,2/7,2/7,4/7 


B13 


t = 32 


1/7,2/7,3/7,3/7 

2/7,2/7,2/7,3/7 


C18 


^(75) 


D23 


j=0 



1/2,1/3,1/3,1/7 29 


J27 


jl4, ^28 


1/3,1/3,1/3,1/3 28 


G36 






H13 


f{-7,7),t^ 


1/3,1/7,1/7,5/7 


18 


i*(0,3,3),t8 


1/3,1/7,2/7,4/7 


B16 


t = 33/2 


1/3,1/7,3/7,3/7 


— 


no covering 


1/3,2/7,2/7,3/7 


— 


no covering 


1/2,1/7,1/7,4/7 27 
1/2,1/7,2/7,3/7 


D46 


i(V-15) 


B2 






125 




1/2,2/7,2/7,2/7 





no covering 


1/3,1/3,1/7,4/7 26 


J12 




1/3,1/3,2/7,3/7 


H51 


j'(4,2), 


2/3,1/7,1/7,3/7 


C12 




2/3,1/7,2/7,2/7 


D47 


t(\/-15) 


1/2,1/3,1/7,3/7 25 


J23 


j" 


1/2,1/3,2/7,2/7 
1/2,1/2,1/7,2/7 24 


H43 


/(-7,0, 14), 


14 


1/3,1/3,1/3,3/7 


D20 


1/3,2/3,1/7,2/7 


B6 


m2(y=3),t=32 




12 


/(0,0,7) 


1/7,1/7,2/7,6/7 


— 


no covering 


1/7,1/7,3/7,5/7 


F24 


i'(\/2T), 


1/7,1/7,4/7,4/7 


G40 




1/7,2/7,2/7,5/7 
1/7,2/7,3/7,4/7 
1/7,3/7,3/7,3/7 

2/7,2/7,2/7,4/7 


— 


no covering 


— 


no covering 


— 


no covering 


— 


no covering 


2/7,2/7,3/7,3/7 


F20 




1/2,1/3,1/3,2/7 23 


J14 


1/2,2/3,1/7,1/7 


J17 




1/2,1/2,1/3.1/7 22 


J25 


7 13, ^26 


1/3,1/3,2/3,1/7 


J15 




1/3,1/7,1/7,6/7 


H36 


i'(9,2),t6 
i = 2"il92/7'i 


1/3,1/7,2/7,5/7 


B34 


1/3,1/7,3/7,4/7 


H52 


J (-3,9) 


1/3,2/7,2/7,4/7 
1/3,2/7,3/7,3/7 


G42 


J (v-7), 




no covering 


1/2,1/3,1/3,1/3 21 


19 


j4(0,3,3), i24 


1/2,1/7,1/7,5/7 


H44 


J^(2,2), 


1/2,1/7,2/7,4/7 


G30 


1/2,1/7,3/7,3/7 


C36 




1/2,2/7,2/7,3/7 


C39 




1/3,1/3.1/7,5/7 20 
1/3,1/3,2/7,4/7 


H42 




All 


t = -l 


1/3,1/3,3/7,3/7 
2/3,1/7,1/7,4/7 


H33 




A14 


t = -1 




B32 


i = 2Vll^ 


2/3,1/7,2/7,3/7 


F23 


j7V2T) 


2/3,2/7,2/7,2/7 


A15 


t = -1 



10 



1/2,1/3,1/7,4/7 19 


Jll 




1/2,1/3,2/7,3/7 


120 


j'(3,7,4) 


^/2,y2}/l,^/7 18 


127 


j5, ^10 


A3 


t = -1 




H37 


i3(9,2),t6 


1/3,1/3,1/3,4/7 


E6 


i"P^(-3, 10) 


1/3,2/3,1/7,3/7 


126 


1/3.2/3.2/7.2/7 
1/7,1/7,1/7,8/7 


— 


no covering 


— 


no covering 


1/7,1/7,3/7,6/7 


C4 


tW2) 


1/7,1/7,4/7,5/7 
1/7,2/7,2/7,6/7 
1/7,2/7,3/7,5/7 


— 


no covering 


— 


no covering 


B31 


t = T^ 


1/7,2/7,4/7,4/7 




no covering 


1/7,3/7,3/7,4/7 


D29 


i(V=3) 


2/7,2/7,2/7,5/7 


— 


no covering 


2/7,2/7,3/7,4/7 
2/7,3/7,3/7,3/7 


— 


no covering 


— 


no covering 


1/2,1/3,1/3,3/7 17 


H49 


j3(-17,51),t6 


1/2,2/3,1/7,2/7 
1/2,1/2,1/3,2/7 16 


133 




122 


f, ^10 


1/3,1/3,2/3,2/7 


A7 


m2(v^),t=-l 




H38 


j3(21,14), t6 


1/3,1/7,2/7,6/7 


G19 


i'(y=3) 


1/3,1/7,3/7,5/7 
1/3,1/7,4/7,4/7 


B17 


t = 3^/2 


— 


no covering 


1/3,2/7,2/7,5/7 


D33 




1/3,2/7,3/7,4/7 


F3 


f{V2) 


1/3,3/7,3/7,3/7 
2/3,2/3,1/7,1/7 


— 


no covering 


A12 


t = -1 




H48 


j3(7,42),t6 


1/2,1/2,1/2,1/7 15 


H45 


i3(-7, 14), 


1/2,1/3,2/3,1/7 


J22 




1/2,1/7,1/7,6/7 


G25 


j2(y^), t4 


1/2,1/7,2/7,5/7 


H40 


j'(2.2) 

t = 3^5/11-'* 


1/2,1/7,3/7,4/7 
1/2,2/7,2/7,4/7 


B33 




no covering 


1/2,2/7,3/7,3/7 


C35 


t{V2l) 


1/2,1/2,1/3,1/3 14 


G39 




H12 






13 


j'(0,0,7),t8 


1/3,1/3,1/3,2/3 


E19 


i''P'(0,28) 


1/3,1/3,1/7,6/7 


D26 


i(V=3) ^ 


1 /'3 1 2 /7 ^ /7 
V-Ji V-Ji / ' : / ' 

1/3,1/3,3/7,4/7 




no covering 




no covering 


2/3,1/7,1/7,5/7 


C27 


t{V7) 

.fiV=7) 


2/3,1/7,2/7,4/7 


G31 


2/3,1/7,3/7,3/7 


C34 


t{V21) 


2/3,2/7,2/7,3/7 


C26 


t{V7) 



1/2,1/3,1/7,5/7 


13 


121 


/(0,4,48) 


1/2,1/3,2/7,4/7 




TTr r\ 

H50 


J (4,2) 


1/2,1/3,3/7,3/7 




C33 


i(\/T3) 


1/2,1/2,1/7,4/7 


12 


G38 


j2(y=7), 


1/2,1/2,2/7,3/7 




F22 


f{V7), 


1/3,1/3,1/3,5/7 




E14 


i"P^(3, 1) 
t = 3^ 


1/3,2/3,1/7,4/7 




B25 


1/3,2/3,2/7,3/7 




B20 


t = 273 


1/2,1/3,1/3,4/7 


11 


H27 


J (-4,4), r 


1/2,2/3,1/7,3/7 




1153 


^•^(-42, 140) 


1/2,2/3,2/7,2/7 




C42 


t(\/385) 


1/2,1/2,1/3,3/7 


10 


H34 


j3(-5,5), t6 


1/3,1/3,2/3,3/7 




D30 


i(^) 


2/3,2/3,1/7,2/7 




C40 


t(\/l05) 


4/3,1/7,1/7,1/7 




E20 


t''P^(21, 14) 


1/9 1/9 1/9 2/7 

/^j /^: /^; / ' 


q 


F,91 


7*'(-24, 34, 141) 


1/2,1/3,2/3,2/7 




115 


1/2,1/2,2/3,1/7 


8 


H47 


j3(7,42),i6 


1/3,2/3,2/3,1/7 




C29 


i(A/7) 
j'(V^), 


1/2,1/2,1/2,1/3 


7 


G35 


1/2,1/3,1/3,2/3 




H35 


j3(0,28), i6 


Table 2.3.8. 








1/8,1/8,1/8,1/8 


36 


F6 
G18 


f{V2), 

-0 / / — . 1 


2/8,2/8,1/8,1/8 


30 


BIO 
G5 


i = 32 

.72(7=1), t4 


1/8,1/8,1/8,3/8 




D24 




1/3,2/8,1/8,1/8 


28 


G8 


j2(y=T), f4 


1/2,1/8,1/8,1/8 


27 


G6 


j'(v^), 


1/3,1/3,1/8,1/8 


26 


C22 


/(-3,2,6), t« 






119 


4/8,2/8,1/8,1/8 


24 


A5 


t = -1 


2/8,2/8,2/8.2/8 

2/8,2/8,1/8,3/8 
1/8,1/8,1/8,5/8 




A18 


i = -1 




— 


no covering 




— 


no covering 


1/8,1/8,3/8,3/8 




B14 
G15 


i = 32 
j'(V^), 


4/8,1/3,1/8,1/8 


22 




no covering 


1/3,2/8,2/8,2/8 




— 


no covering 


1/3,2/8,1/8,3/8 

1/2,2/8,2/8,1/8 




116 


i*(2,8,8) 


21 




no covering 


1/2,1/8,1/8,3/8 




G16 


j2(y=2), f4 


1/3,1/3,2/8,2/8 


20 


A9 
H22 


t = -1 
j3(5, 10), i6 


1/3,1/3,1/8,3/8 
2/3,2/8,1/8,1/8 




D16 


i(v^) 




A13 
G3 


m2(v^),i=-l 

i'(y=T), 


1/2,1/3,2/8,1/8 


19 


J21 





11 



1/2,1/2,1/8,1/8 


18 


F4 


j2(V2), i4 

/,tl2 






Jl 


4/8,4/8,1/8,78 




— 


no covering 


4/8,2/8,2/8,2/8 






no covering 


4/8,2/8,1/8,3/8 




Bl 


t = 22 


1/3,1/3,1/3,2/8 




E5 


f''P'(-3, 10) 
j2(y=2), f4 


1/3,2/3,1/8,1/8 




G12 


2/8,2/8,1/8,5/8 




Dl 


t(x/=T) 


2/8,2/8,3/8,3/8 




C5 


t(\/3) 


2/8,6/8,1/8,1/8 




— 


no covering 


1/8,1/8,1/8,7/8 
1/8,1/8,3/8,5/8 




— 


no covering 






no covering 


1/8,3/8,3/8,3/8 




— 


no covering 


1/2,1/3,1/3,1/8 


17 


J20 


j9,tl8 


4/8,1/3,2/8,2/8 
4/8,1/3,1/8,3/8 


16 


— 


no covering 




— 


no covering 


1/3,1/3,1/3,1/3 




D8 
F2 




1/3,2/8,1/8,5/8 
1/3,2/8,3/8,3/8 




B30 


t = 22132 




~ 


no covering 


1/3,6/8,1/8,1/8 




G22 




1/2,4/8,2/8,1/8 


15 


Gl 


i = 34 


1/2,2/8,2/8,3/8 




B15 


1/2,1/8,1/8,5/8 




G9 


j2(y=T), t4 


1/2,1/8,3/8,3/8 




C8 


j^(x/2), t4 


4/8,1/3,1/3,2/8 


14 


F8 


4/8,2/3,1/8,1/8 




— 


no covering 


1/3,1/3,1/8,5/8 




H23 


j3(5,10),t6 


1/3,1/3,3/8,3/8 




C28 




2/3,2/8,2/8,2/8 






no covering 


2/3,2/8,1/8,3/8 




B29 


t = 5'^/22 


1 /O /I /O 1 /O 1 /o 

1/2,4/8,1/3,1/8 


13 


T1 O 

118 


J (-3,2,6) 


1/0 i/Q 2/a 3/a 




xizy 


/(-1, 2) 


1 /O 1 /O 2 /C 2 /Q 

-^/z,YZ,^/o,''/o 


1 

Iz 


v^l 


/o^ 


I/O I/O 1/C "^/Q 




jjy 
G14 


/ o^ 
Hv -^) 

J (V-2), t 


4/8,1/3,1/3,1/3 


12 


E2 


t''Pi(3,2) 


1/3,2/3,2/8,2/8 




C21 


t(V6) 
J^(a/=2) 


1/3,2/3,1/8,3/8 




GIO 


1/2,1/3,1/3,3/8 


11 


F25 


j2(V22), f4 


1/2,2/3,2/8,1/8 

1 /O 1 /O 1 /o o /o 

V2,V2,V3,2/8 




110 


j4(-2,4,-l) 


10 


H20 


/(5,10),i6 

■2/ /n\ J.4 

.7 (v6), 


1/3,1/3,2/3,2/8 




F19 


o/o O/O 1 /O 1 /o 

2/3,2/3,1/8,1/8 




G32 


1 / /-t f\\ 
HV 10) 

i2(^), fl2 


1/2,1/2,1/2,1/8 


9 


G17 


1/2,1/3,2/3,1/8 




J2 




1/2,1/2,1/3,1/3 


8 


F5 


f{V2), ^4 
i'(^/=2), 






G13 


1/3,1/3,1/3,2/3 




D13 


t(x/=2) 



Table 2.3.9. 

3/9,1/9,1/9,1/9 

1/9,1/9,2/9,2/9 

1/3,1/9,1/9,2/9 
1/2,1/9,1/9,1/9 
1/3,1/3,1/9,1/9 

3/9,3/9,1/9,2/9 
3/9,1/9,1/9,4/9 

3/9,2/9,2/9,2/9 

6/9,1/9,1/9,1/9 
1/9,1/9,2/9,5/9 
1/9,2/9,2/9,4/9 
1/3,3/9,3/9,1/9 
1/3,3/9,2/9,2/9 
1/3,1/9,1/9,5/9 
1/3,1/9,2/9,4/9 
1/2,3/9,1/9,2/9 
1/2.1/9,1/9.4/9 
1/2,2/9,2/9,2/9 
1/3,1/3,3/9,2/9 
1/3,1/3,1/9,4/9 
3/9,2/3,1/9,1/9 
2/3,1/9,2/9,2/9 
1/2,1/3,3/9,1/9 
1/2,1/3,2/9,2/9 
1/2,1/2,1/9,2/9 
1/3,1/3,1/3,3/9 
1/3,2/3,1/9,2/9 

1/2,1/3,1/3,2/9 
1/2,2/3,1/9,1/9 
1/2,1/2,1/3,1/9 
1/3,1/3,2/3,1/9 
1/2,1/3,1/3,1/3 

Table 2.3.10. 



24 — 
A24 
H8 
22 H9 
21 G52 
20 HI 
J3 
18 — 



E3 
C23 

16 A6 
07 
H25 
H6 

15 H4 
D49 

14 D39 

G21 

D2 

09 

13 117 
C41 

12 H7 
D19 
B7 
H2 

11 H32 
H31 

10 J4 
H3 
9 G26 



no covering 
t = -1 
j'(-3,4),t6 
/(-3,4),t6 
j2(y=T5), ti2 

i'(-3,i) 



.12 



no covering 
no covering 
no covering 

i'^P'(3,2) 
t{V6) 

no covering 
t = -1 
t(\/3) 

J^(3,2) 
i'(0,2) 
i(y=39) 
no covering 

7) 



3),i4 



fi 
j"(0,6,9) 



i(V273) 

j = 
t = 3^ 

.7^(0,3) 
/(6, 1), 

i2(^),fl2 



1/10,1/10,1/10,1/10 



5/10,1/10,1/10,1/10 

2/10,2/10,2/10,2/10 
2/10,2/10,1/10,3/10 
2/10,4/10,1/10,1/10 
1/10,1/10,3/10,3/10 
1/3,2/10,2/10,2/10 
1/3,2/10,1/10,3/10 
1/3,4/10,1/10,1/10 
1/2,2/10,2/10,1/10 
1/2,1/10,1/10,3/10 



24 D6 
F17 
18 E7 

B24 

D4 

025 

16 — 
02 
D42 

15 015 
D44 
D50 



t{V-i) 

f{V5), i4 
t'P^S, 10) 
no covering 
t = 273 
^(^/^) 

no covering 

t{V-u) 

t(V-15) 
t(\/=39) 



12 



1/3,1/3,2/10,2/10 

1/3,1/3,1/10,3/10 

2/3,2/lO,l/lO,l/lO 

1/2,1/3,2/10,1/10 

1/3,1/3,1/3,2/10 
1/3,2/3,1/10,1/10 
1/2,1/3,1/3,1/10 
1/3,1/3,1/3,1/3 

Table 2.3.11. 



14 C38 
D48 
F14 
13 J5 
12 FIO 
G7 
E13 
H30 
11 J8 
10 B19 



t(V21) 
t"Pi(3,l) 

j3(0,10),i6 



25/5 



1/11,1/11,1/11,4/11 
1/11,1/11,2/11,3/11 

1/11,2/11,2/11,2/11 
1/3,1/11,1/11,3/11 
1/3,1/11,2/11,2/11 
1/2,1/11,1/11,2/11 
1/3,1/3,1/11,2/11 
2/3,1/11,1/11,1/11 
1/2,1/3,1/11,1/11 
1/3,1/3,1/3,1/11 

Table 2.3.12. 

3/12,1/12,1/12,1/12 
2/12,2/12,1/12,1/12 
1/3,2/12,1/12,1/12 

1/2,1/12,1/12,1/12 
1/3,1/3,1/12,1/12 

Table 2.3.13. 

1/13,1/13,1/13,2/13 
1/3,1/13,1/13,1/13 

Table 2.3.14. 
1/14, 1/14, 1/14, 1/14 

Table 2.4.5. 
1/5,1/5,1/5,1/5 



18 



D15 



16 D41 

C13 
15 H26 
14 H28 

E15 
13 J9 
12 G45 



no covering 

no covering 
t{V=7) 

t{V3) 

j3(-4,4),t6 
i'(-4,4),t6 
t"P'(-ll,22) 

^12 



1/2,1/4,1/5,1/5 
2/4,1/4,1/4,1/5 
1/5,1/5,1/5,3/5 
1/5,1/5,2/5,2/5 
1/2,1/4,1/4,1/4 
2/4,1/5,1/5,2/5 
1/4,1/4,1/5,3/5 
1/4,1/4,2/5,2/5 
1/2,1/4,1/5,2/5 
1/2,1/2,1/5,1/5 

2/4,2/4,1/5,1/5 

2/4,1/4,1/4,2/5 
1/4,3/4,1/5,1/5 
1/5,1/5,1/5,4/5 
1/5,1/5,2/5,3/5 
1/5,2/5,2/5,2/5 
1/2,2/4,1/4,1/5 
1/2,1/2,1/4,1/4 



18 D21 j = 

Bll i = 32 

16 D14 t(v^) 

G24 f{y^),t^ 

15 G27 j2(y33)^ii2 

14 F9 fiV^) 
G23 

18 E16 t"Pi(-l,2) 

16 G28 j2(y=3),ii2 



18 C2 t{V2) 



2/4,2/4,1/4,1/4 
2/4,1/5,1/5,3/5 
2/4,1/5,2/5,2/5 
1/4,1/4,1/4,3/4 
1/4,1/4,1/5,4/5 
1/4,1/4,2/5,3/5 
1/2,1/4,1/5,3/5 
1/2,1/4,2/5,2/5 
1/2,1/2,1/5,2/5 
2/4,2/4,1/5,2/5 

2/4,1/4,1/4,3/5 
1/4,3/4,1/5,2/5 
1/2,2/4,1/4,2/5 
1/2,3/4,1/5,1/5 
1/2,1/2,2/4,1/5 
1/2,1/2,1/2,1/4 

Table 2.4.6. 



1/4,1/4,1/5,1/5 
1/4,1/4,1/4,1/4 



1/5,1/5,1/5,2/5 
2/4,1/5,1/5,1/5 
1/4,1/4,1/5,2/5 



24 A20 
F12 
G48 
22 130 
20 C17 
D45 
Ell 
E9 



j2(y=45), t4 
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G51 



t(V-15) 
i"P'(5,10) 
rP\5, 10) 
no covering 



17 J18 
16 H18 
E18 
H16 
15 E25 
14 D38 
CIO 
H41 
13 131 
12 F16 
111 
A2 
G50 
A4 
F15 
D3 
B23 

11 129 

10 C14 
C30 
H17 
F13 
D43 
C16 

B8 

9 124 
D27 
8 H19 

C19 
B26 
7 H39 
G29 
6 G46 
5 ElO 



f, tie 

j3(5,10), i6 

i"P'(0.10) 
f{5, 10) 

t{^) 
f 

t = -1 
j2(y=T5), 

m?{^/^),t=-l 

t = 273 

no covering 

f 

tiVio) 

j3(5, 10), 

i(A/=T5) 

no covering 

t = 3^ 

no covering 

f 

j3(5,10), i6 
no covering 

t{V6) 
t = 372^ 
j'(2,2) 

t''P^(5,10) 



1/6,1/6,1/6,1/6 
1/4,1/4,1/6,1/6 

3/6,1/6,1/6,1/6 
2/6,2/6,1/6,1/6 
1/4,1/4,1/4,1/4 



2/4,2/6,1/6,1/6 
3/6,1/4,1/4,1/6 
2/6,2/6,1/4,1/4 



16 B28 
Dll 

14 Cll 
113 

12 D18 
B5 
B4 
B22 



10 



H5 
D32 



t = 3^/2^ 

i4(0,4,12),t8 

J=0 

t = 32 

t = 22 

t = 273 

no covering 

j'(O^), 
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k, £, Tfi 


Max. 


Br. patterns 


Coverings 


Moduli field dei 


^ree 










degree 


Total 


No GOV. 


Orb. 


Total 


1 


2 


3 


4 


5 


6 


> 7 


2, 3, 7 


60 


152 


30 


140 


427 


51 


27 


25 


11 


8 


5 


13 


2, 3, 8 

■"7 '-'7 ^ 


36 


65 


16 


58 


130 


23 


24 


4 


4 




2 


2 


2, 3, 9 


24 


32 


6 


29 


67 


12 


3 


11 


1 




2 




2, 3, 10 


24 


20 


2 


21 


38 


13 


5 


1 






2 




2, 3, 11..14 


18 


18 


2 


18 


33 


9 


6 


2 






1 




2, 4, 5 


24 


40 


5 


42 


91 


21 


10 


6 


1 


4 




1 


2, 4, 6/7/8 


16 


24 


5 


25 


43 


16 


4 


1 


4 








2, 5/6/7 


12 


12 


1 


15 


22 


10 


3 


2 










3, 3/4/5 


12 


15 


2 


16 


21 


12 


3 


1 










Total 




378 


69 


366 


872 


167 


85 


53 


21 


12 


12 


16 



Table 1: Statistics of Belyi maps 



Table 2.4.6 (continued) 



Table 2.5.6. 



1/2,2/6,1/4,1/6 9 II i4(0,8,12) 
1/2,1/2,1/6,1/6 8 C3 t{V2) 

Gil fiV=^), 

2/4,2/4,1/6,1/6 — no covering 

2/4,2/6,1/4,1/4 D12 t(v^) 

1/4,3/4,1/6,1/6 G20 f lV^), 

1/2,2/4,1/4,1/6 7 112 j4(0,4,12) 

1/2,1/2,1/4,1/4 6 C24 f(V6) 

D25 t{^/^) 

Table 2.4.7. 

1/7,1/7,1/7,2/7 12 — no covering 

2/4,1/7,1/7,1/7 10 — no covering 

1/4,1/4,1/7,2/7 D34 

1/2,1/4,1/7,1/7 9 17 j4(2,8,9),t8 

2/4,1/4,1/4,1/7 8 G37 

1/2,1/4,1/4,1/4 7 E24 rP\-7,U) 

Table 2.4.8. 

1/8,1/8,1/8,1/8 12 — no covering 
1/4,1/4,1/8,1/8 10 B9 t = 32 

G4 .f{V^),t^ 
1/4,1/4,1/4,1/4 8 A16 t = -1 

Table 2.5.5. 

1/5,1/5,1/5,1/5 12 A19 t = -1 

Fll f{V5) 

G47 j2(y=T5), ^4 

1/5,1/5,1/5,2/5 10 E8 t''Pi(5,10) 

1/5,1/5,1/5,3/5 8 E17 t'^P'^O, 10) 

1/5,1/5,2/5,2/5 H15 j3(5_J_0) 

1/2,1/5,1/5,2/5 7 D37 t{V^) 

1/2,1/2,1/5,1/5 6 Al t=-l 

G49 f{y/^), 



1/6,1/6,1/6,1/6 10 C31 t(VlO) 

2/6,1/5,1/5,1/5 8 — no covering 

1/5,2/5,1/6,1/6 D36 t{V^) 

1/2,1/5,1/5,1/6 7 H24 f{3,l),t^ 

Table 2.5.7. 



1/5,1/5,1/5,1/7 J 


i E23 t^P'(2,2) 


Table 2.6.6. 




1/6,1/6,1/6,1/6 J 


? B27 t = 3^/2^ 




DIO i(\/^) 


Table 3.3.4. 





1/4,1/4,1/4,1/4 


12 


A17 


t = -1 


1/3,1/3,1/4,1/4 


10 


A8 

H21 


t = -l 

,73(5, 10), i6 


2/4,1/4,1/4,1/4 


9 




no covering 


1/3,1/3,1/3,1/4 




E4 


t"P'(-3,10) 


1/3,1/3,1/3,1/3 


8 


D7 


^(^/=2) 






Fl 


m*, i(\/2),i4 


2/4,1/3,1/3,1/4 


7 


F7 


2/4,1/3,1/3,1/3 


6 


El 


t"PH3,2) 


1/3,2/3,1/4,1/4 




C20 


i(V6) 
i'(^/6), 


1/3,1/3,2/3,1/4 


5 


F18 



Table 3.3.5. 



1/5,1/5,1/5,1/5 
1/3,1/3,1/5,1/5 
1/3,1/3,1/3,1/5 
1/3,1/3,1/3,1/3 

Table 3.4.4. 



9 — no covering 

7 C37 t{^/2i) 

6 E12 f^Pi(3,l) 

5 B18 i = 2V5 



1/4,1/4,1/4,1/4 
1/3,1/3,1/4,1/4 



6 B3 i = 22 

B21 t = 2^/3 
5 D31 i(\/^) 
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• if j £ Q \ {0} and t G Q, a value of t is given in the factorized form (as 
motivated by j|E|; 

• if j G Q \ {0} and the t-field is quadratic, t{^/a) means the field Q{\/a); 

• if J e Q \ {0} and the t-field degree is greater than 2, t^P^{a, b) means the 
splitting field (of degree 6 in our examples) of the polynomial + aX + b. 

The size of the Galois orbil[^is equal to the degree of the moduli field. In most 
cases, the j/-ficld and the moduli field coincide]^ The moduli field is indicated 
only if it differs from the j-field, either with m?{^/a) if its degree is 2, or with 
to'' if its degree d> 2. 

More information about each computed Galois orbit can be found in the 
tables of Appendix ^ and our website [17]. Statistics of computed Belyi 
coverings is presented in Table [T] There are no coverings for 69 out of 378 
derived branching patterns. 



4 The moduli field and obstruction conies 

In this section, let O denote the group of Mobius transformations: 

( ax + h 1 

0= I |a,6,c,de Q with ad-bc^Q \ = Aut(Q(a;)/Q). 

I C-X I d J 

Two rational functions Lpi,ip2 G Q(a;) are called Mdhius-equivalen^ denoted 
(pi ~ (p2j if there exists ^ G O with (pi o fi — ip2- 

A realization field of a Belyi covering tp is a number field over which a Mobius 



equivalent function o /i is defined. The r-field from Definition 2.3 is such a 
field, but often not of minimal degree. This section describes how, for each Belyi 
function ip in our table, we determined its moduli field and realization fields. 

Definition 4.1. Let ip £ Q{x) he a Belyi function. Let G = Gal(Q/Q). The 
moduli field M^p is the fixed field of {a G G\p ^ a{p)}. 

Clearly C K^p for any explicit p G K^{x) over some realization field 
Kip. The moduli field is known to be equal the intersection of the realization 
fields of p. It can be computed directly from the given definition, by checking 
which Galois conjugates of p> are Mobius-equivalent to p. The computed Belyi 
functions p always had : Q(j)] < 2, where Q(j) is the j-field. But Q(j) C 
C K^. Therefore = if = Q(j). For those ^ Q(j), let a be 
the non-trivial element of G&\{K^/'^{j)). Our moduli fields are then determined 
by just checking whether p ^ <y{p)- 



^Our notation allows to count the total number of dessins d'enfant in selected Galois 
orbits rather quickly in Tables 2.3.7-3.4.4. Each fourth column starts either with the m'* or 
j'^ notation (where d is the size of the Galois orbit), or starts with a symbol without a numeric 
power (in which case d = 1). 

* The moduli field always contains the j'-field, as the j'-value is an invariant of Mobius 
transformations. In our encountered cases, the moduli field is at most a quadratic extension 
of the j-field. 

^To determine whether ipi ^ ip2, we can factor the numerator of fiix) — ip2{y). If it has a 
factor p{x,y) S Q[a:,3/] with deg^(p) = dcgj^(p) = 1, then ipi ~ ip2 and one finds fi by solving 
p{x, y) = with respect to y. 
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Id 


branching 


d 


[fc^rn] 


Moduli 


Obstruction 


Bad 




fractions 






field 


conic 




primes 


B12 




36 


[237] 


Q 


2 

u 


2 

+ W + 


7 


7, 00 


Co 


73, 73, 2/7, 2/7 


32 


[237] 


Q 




+ + 


1 
1 


2, 00 




1/0 1/0 1 //I 1 //f 


1 n 


[^45J 


(TTl 
V 


„,2 
U 


0„,2 

+ zw - 




5, (X) 


T\A ^ 
U40 


1 //I 1 //I 1 //I 1 //I 

74, 74, 74, 74 


on 


[^45J 


(TTl 
V 


„,2 


0„,2 

+ IV - 




5, 00 


i 


1/01 /o 1/0 1/0 

73, 73, 73, 73 



» 


[334] 


<lj(v2) 




+ Sw'^ - 


-f V z — i 


3, 00 


F4 


1/0 1/0 1/0 1/0 

72,72,78,78 


18 


[238] 


Q(V2) 




9 


-1 

1 


00, 00 


F6 


1/8, 1/8, 1/8, V8 


36 


[238] 


Q(y2) 


9 


+ + 


1 


ex, 00 


Fll 


1/5, 1/5, V5, V5 


12 


[255] 


Q(%/5) 




+ 21-2- 


f V5 


5, 00 


HI 


1/3, 1/3, 1/9, 1/9 


20 


[239] 


Q(ReC9) 


U2 




ReCg 


00, 00 


HIO 


1/2, 1/2, 1/7, 1/7 


30 


[237] 


Q(ReC7) 


U2 


+ - 


ReC7 


cx, cx 


Hll 


1/3, 1/3, 1/7, 1/7 


44 


[237] 


Q(ReC7) 


U2 


+ _ 


ReC7 


00, 00 


H12 


1/2, 1/2, 1/3, 1/3 


14 


[237] 


Q(ReC7) 


U2 


+ - 


ReC7 


cx, 00 


H13 


1/3, 1/3, 1/3, 1/3 


28 


[237] 


Q(ReC7) 


u2 


+ _ 


ReC7 


00, cx 


H14 


1/7, 1/7, 1/7, 1/7 


60 


[237] 


Q(ReC7) 


u2 


+ - 


RcC7 


cx, cx 



Table 2: Belyi functions with an obstruction. 



Among the minus-4-hyperbolic Belyi functions, there are 14 Galois orbits 
for which the moduli field not a realization field. They are given in Table [2j 
The column "Obstruction conic" gives an answer to the question: What are the 
realization fields of a Belyi map Lpl This question is illuminated in j ]4.1| and j ]4.2[ 
The last column gives an alternative characterization of obstruction conies, as 



explained in [4.3 



Two Belyi functions are Mobius-equivalent if and only if they have the same 
dessin d'enfant up to homotopy. Thus, the moduli field of a dessin d'enfant is 
well defined. The number of different dessins (up to homotopy) in a Galois orbit 
is equal to the degree of the number field. The dessins d'enfant of most of the 
Belyi maps of Table [2] are depicted in Figure [3j The other Galois orbits with 
obstructed dessins can be found in Figures [T] and [4j 

The interesting questions whether a dessin has a moduli field C M, and if 
so, does it have a realization over M, are considered in [5]. Although all moduli 
fields in the obstructed cases are real, not all their dessins have a reflection 
symmetry (i.e., have a realization over M). Rather, their complex conjugates 
are equivalent to the original up to homotopy that permutes the cells, reflecting 
a non-trivial Mobius equivalence. The number of these skew- symmetric dessins 
depends on the number of bad c»-primes shown in the last column of Table [2j 
The moduli field for HI, HIO - H14 has three infinite primes, but only two of 
them are bad. Therefore one dessin in those orbits has a reflection symmetry, 
and the other two are skew symmetric. Likewise, Fl and Fll each have one 
dessin with M-realization and one dessin without. 

4.1 Obstructions on realization fields 

If the moduli field M^p is a realization field, the other realization fields are 
extensions of M^. Otherwise, the realization fields are determined by a conic 
obstruction. For each of the cases of Table [2j the realization fields are those 
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extensions of that have a rational point on the conic curves given in the 
sixth column. 

If </3 e Q(x), we denote 



0^ = {^l(^0\^o^, = ^}^ Aut(Q(a:)/Q((^)). 
We have \0^\ > 1 in the following four cases: D45, F6, H13, H14. Then 



\0^\ = 2, and the dessins have a non-trivial automorphism. 

Lemma 4.2. If ip £ K{x) is a Belyi function, and |0y| — 2, then we can write 
if = f o f2 for some /, /2 G K(x) with f2 of degree 2. If ip has a realization over 
L, then so does f. 

Proof: Let n be the non-identity element of O,^. Since ip is invariant under 
Ga.\{K/K), the same must be true for n, and so G K{x). Since n has order 
two in Aut(Q(a;)/Q), its fixed field F C K{x) has index 2. Define /2 € F as x/x 
if Xfx ^ K, and x + jj, otherwise. Then /2 has degree 2, and hence it generates 
F. Since </9 e F, it follows that ip = f o f2 iov some / G i4r(a;). 

If Lp has a realization G -^^(a;), then from the non-identity element of 
we can compute explicit /2 and / in L{x) in exactly the same way. So if ip> has 
a realization over L, then so does /. □ 

Lemma 4.3. Let L he a field of characteristic 0. Then 

C30 has a realization over L D45 has a realization over L. 

Proof: "=>": Let pi S i(a;) be a realization for C30. The branching structure 
of C30 above 1 is 1,1,2,2,2,2 and so (^"^({1}) C QlJ{oo} contains 6 points, 
two of which (call them Pi,P2) are unramified. Let S := {Pi,P2}. Since 
ip is invariant under Gal(L/L), the same must be true for S. So either 
Pi,P2 G L lJ{oo}, or Pi, P2 are the roots of an irreducible quadratic polynomial 
Q S L[x]. In both cases one can construct an explicit /2 G L{x) of degree 2 that 
branches above Pi and P2- Then ipo f2 G L{x) is a Belyi function with the same 
branching type as D45. The Heun table shows that among the Belyi functions 
with this branching type, the only one that allows C30 as a decomposition factor 
is D45. Hence (p o f2 is a realization of D45 over L. The converse statement 
follows from Lemma l42l □ 

The same proof also works for F6, F4, for H13, H12, and for H14, HIO. So 
it suffices to answer Question (3) for the cases with \0,p\ = 1. We shall do this 
for an example, namely C30. The same works for the remaining cases, and the 
results are listed in Table [21 

4.2 The obstruction conic for C30 

The moduli field for C30 is Af = Q. Our website [17j lists this Belyi function 

_ 2((a;2 + 5)^/^ - 3x^ - 60x + 15){x^ + 5x - 5)^ 



in K{x) where K = Q{^/—3). Suppose ^ ^ g for some g e L{x) with \/— 3 ^ L. 
Write ip ^ g o fi for some /i G O, which must be unique because Oip is trivial. 



(p = 



(12a;)5 
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That implies ^ £ L(-\/— 3)(a;) since (^,17 G 3)(x). So we can write /i = 

(ax + &)/ (ca; + d) with a, b,c,d £ L{\/—3) and ad — be 0. 

Without loss of generality we may assume c = or c = 1. The case c = 
is a quick computation, and does not lead to solutions (with ad — be ^ and 
^ L) so we set c := 1. Write a = og + ai\/^, 6 = 6q + 61 -^/^ and 
d — d(j + 3, for some oq, ai, fopi ^ij c^Oi "^i ^ We can replace /i by ^ — ag, 

this corresponds to replacing .g(a;) by (/(x + ao), which is still in L[x). After this, 
we may assume ag = 0. The case ai = does not lead to solutions so we may 
suppose oi 7^ 0. Replacing by ^/ai corresponds to replacing g{x) by g{aix), 
after this, we may assume ai = 1. Then we can write ip = g o fi where 



+ 60 + &i\/-3 
a; + do + di\/— 3 



(4.1) 



for some bo,bi,do,di € L. 

Let cr : "v^— 3 1— > — -^/^ be the non-trivial element of Gal{L{^/~3) /L). Now 
a{(p) ~ (/9 (otherwise the moduli field would have been K). We compute the 
unique v G O with a{(p) = ipo u hy factoring the numerator of ip{y) — criip{x)), 
and find 1/ = — 5/x (uniqueness implies v € K{x). Now 

9 ° f^(Ai) = cr(5) ° cr(^) = a{g o fi) = cj{ip) = ^pov = go fiov 



and hence (t(/i) = o j/ (using again \0g 



1). Write the numerator 



of cr(^) - fiov a.s X^Lo I]i=o C'jj a;*(V^)^ with dj € Q[6o, do, ^i]- The 
equation (T(/i) = fio v corresponds to setting all Ci_j to 0, resulting in: 

^1 + do = 60 ^ 3di = 5 + ^oc^i ^ ^ic^o = 



which reduces to 



bi 



-do, 60 = 3di, and dp + 3d]^ + 5 = 0. 



(4.2) 



So if has a realization over L then dQ + 3d^ + 5 = has a solution do, di G L. 
Conversely, for any solution do,di e L, one obtains an explicit realization of 



C30 over L, namely g — (p o ^, with /i as in (4.11,(4.2) 



4.3 Cohomological obstructions 



The conic (4.2 ) describes all fields L over which C30 has a realization. The conic 
+ 2t^ + 5 = is birationap^ to ( |4.2[ ) and hence it describes the same fields 



C30 has an obstruction; it has no realization over its moduli field Q because (4.2 1 
has no Q-rational points. Table [2] lists all entries of the Heun table that have 
an obstruction. 

Obstructions can be described in several ways: skew fields, 2-cocycles, or 
bad primes. If p is a prime of the number field M^, then p is called a bad 
prime if if has no realization over the completion of at p (M is denoted by 
p = 00). The set of bad primes completely describes the obstruction (a conic is 
determined up to birational equivalence by its bad primes). 



^''We have an implementation 17 to find birational maps for conies over Q. For number 
fields, we followed the same method, computing F1,F4,. . .,H14 case by case. 
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li If G Q{x) has \0,p\ = 1 then (p defines a 2-cocy cle in a natural way, 
representing an element of 7?^(Gal(Q/Q), O). Section 4.2 spells this out in 



an explicit example; we computed the image of (id, a) under the 2-cocycle (the 
v G O for which id((p) o = a{(f)), and then used that to determine the fields 
L for which this 2-cocycle maps to the identity in i?^(Gal(Q/L), O). 

A conic over with a point defined over an odd-degree extension of Af^ 
will necessarily have a point over M^. So an obstruction can only occur for 
branching patterns where in each branching index appears an even number of 
times in each partition. In particular, the branching fractions in Table [2] must 
be of the form a, a, b, b for some a, b. 

While the conies for C6 and F4 look the same, their sets of bad primes differs. 
The explanation is that the moduli field of F4 has two infinite primes (two real 
embeddings) while the moduli field of C6 has only one real embedding. 

4.4 Conic-realizations 

We say that a Belyi function ip has a conic-realization if there exist a conic C 
defined over M^p such that the function field of C over Af^ has a Belyi function 
f : C ^ with the same dessin as (p. 

In the example C30, let L = Q(do)Mi]/(rf^ + 3d? + 5) be the function field 



of the conic. Then (p ^ g where g := (p o fjr^ g L{x) with /i as in (4.1 1,(4.2). 
Substituting an arbitrary rational number for x in g produces an element of L. 
This element turns out to be a conic- realization for C30. In the same way we 
can find a conic-realization for all cases with jO^j — 1. 

Given a realization for C30 over L, we can use lemma |4.3| to obtain a 
realization h € L{x) for D45. However, if we substitute a rational number 
for X, the resulting element of L will not have the D45 dessin, instead, it has the 
C30 dessin. It follows from jlj that the |0i^| = 2 cases can not have a conic- 
realization (in fact, [5, Section 4] explains why no element of M(7i)[w]/ [v? + 1) 
can have D45 or F6 as a dessin). 

Sometimes a conic-realization is a compact way to represent a Belyi function, 
for example: 



where Ei = + v"^ + 2, 



4 

C30: («-«-10) ( ^ + ^) eQ(u)M/(u2 + 3t'2 + 5), 



1 ^ 5 
li — 1 



and 

Fll : v\ ) e Q{V5){u)[v]/{E2) 

where E2 = 29u'^ -Auv + v'^ ~Au + 1102v + l + 2{6u'^ ~uv-u + 2A6v)\/5 which 
is much shorter than any (p e Q(a;) we found for Fll. 
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Fl 




F7 




Fll 





or 





Figure 4: Indefinite cases of moduli fields 



4.5 Ambiguous moduli fields 

The moduli field for the Galois orbit F7 is M = Q(\/3 + 6^2) by the standard 
definitions. However, the branching pattern is [4] + 2 + 1 = 2 [3] + 1 = 2 [3] + 1 
has two symmetric fibers 2 [3] + 1. By using a Mobius transformation on Pj, a 
Belyi function for F7 can be expressed over Q(\/2): 

(x-\/2)(l62;^-32x^ + (26-7\/2)x-14+10V2)(l6a;^(a;+l-V2) + (58-41\/2)(2;-l)) 

The branching fibers for this rational function are z ~ oo and 



(4.3) 



753 - 531V2 
128 ' 



(4.4) 



with two fibers conjugated in M D Q(a/2). Defining a Belyi function 
by requiring branching in (at most) any 3 fibers, not specifically {0,1, oo}, 
would make no geometrical difference because of Mobius transformations. But 
evidently, there are arithmetic consequences for moduli and realization fields. 

The number of dessins for F7 is 2 or 4 depending of whether the dessins 
are counted up to Mobius equivalence on or not. Figure |4] depicts two of 
the dessins for F7. The other two are obtained by swapping the color labeling 
of black and white vertices, which corresponds to swapping the two symmetric 
fibers. If the symmetric fibers are put at z = 0, z = 1, the transformation 
z I— ?> 1 — z swaps the two symmetric fibers and changes the sign of y/i + 6%/2. 
The second depicted dessin for F7 is realized over K. 

Even more interestingly, the Galois orbits Fl and Fll demonstrate an 
interesting mix of a conic obstruction and ambiguous moduli field. Their 
realization fields are obstructed by the conies in Table[2]if we insist in having the 
branching fibers at {0, l,oo}. But the Belyi functions can be expressed (after a 
Mobius transformation on V\) over their moduli fields if the location of fibers is 
not fixed. In particular, the following rational function in Q(\/2)(x) gives Fl: 

3:° + (10-272)z-'' + (2-3572)a'^ + (176+8472)3:^ + (66+6872)a:^+(44+2872)a:-20-1472 
(2;2+2+\/2) (2;2-2\/2x-2-\/2)'' 



The ramified fibers are z = oo and z = ±| Reflecting this, the first 

dessin of Fl in Figure |4] is symmetric if vertex coloring is ignored, but the black 
and white vertices are interchanged by the complex conjugation. The dessins 
for Fll are drown in two variations: first compactly, by hiding white vertices of 
order 2; then assigning the black and white vertices to represent points of order 
5 to show the fiber interchanging symmetry. 
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Examples of coverings with this dual interpretation of the moduli field are 
given in ^ . The examples of Pharamond have the branching patterns 4 + 2 = 
4+1 + 1 = 4 + 1 + 1 or 4 + 2 + 1 = 4 + 2 + 1 = 4 + 2 + 1. The moduli field 
for the first pattern is Q(-\/5), but the rational function can be written over Q 
if the fiber location is not fixed. There are two orbits for the pattern with three 
symmetric fibers. A rational function for the first orbit can be similarly written 
in Q(a/2), though the moduli field is Q{\/—1 — 2\/2) by definition. A function 
for the other orbit can be written in Q(-\/— 6), while the moduli field is of degree 
12, obtained by adjoining the roots of the polynomial — z'^ + {3 + ^/^)z — 3. 



5 Application to Heun functions 

The minus-4-hyperbolic Belyi functions have application to transformations 
between hypergeometric and Heun functions (or their differential equations). 
This allows to express some Heun functions in terms of better understood 
hypergeometric functions. In fact, we utilize this application in our algorithms 
to compute the Belyi functions. 

The Gauss hypergeometric equation 



dz^ 



A + B-C+l\ dy{z 



z-l 



AB 



dz 



z{z-l) 



v{z)=Q 



and the Heun differential equation 
d'^y{x) , f c d a + b-c-d+l 



dx"^ 



X- 1 



X — t 



\ dy{x) abx — q 

J dx x{x — l){x — t) 



(5.1) 



V{x) = 



(5.2) 

are second order Fuchsian equations [21] with 3 or 4 singularities, respectively. 
The singular points are z = 0, 1, oo and a; = 0, 1, oo. If C ^ Z, a basis of local 
solutions of (5.1) at a; = is given by the famous Gauss hypergeometric series: 

(5.3) 



2F1 



A, B 
C 



2F1 



l + A-C,l + B-C 
2-C 



The starting powers 0, 1 — C of the local parameter z are the local exponents at 
z = Q. The local exponents at z = 1 are Q,G — A— B, while the exponents at 
z = 00 are A, B. The local exponents for Heun equation (5.2) are 



at a; = : 0, 1 — c; 
at a; = 1 : 0, 1 — d; 



at a; = cx) : a, 6; 
ai X = t : Q,c + d ~ a — b. 



The local solution at a; = with the exponent is denoted by 

Hn( 





a, b 




c 








c; d 





(5.4) 



The parameter q is an accessory parameter^ it does not influence the local 
exponents. If c ^ Z, then an independent local solution at a; = is 



Hn 



a — c+1, 5 — c 
2 - c; d 



1 



(5.5) 
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with = g - (c - l)(a + b- c- d + dt+1). 

A pull-back transformation of a differential equation for y{z) in d/dz has the 
form 

z I — > Lp{x), y{z) I — > Y{x) = 6'(x) y{ip{x)), (5.6) 

where v^C^^) is a rational function, and ^(a::) is a radical function (an algebraic 
root of a rational function). Geometrically, the transformation pulls-back a 
differential equation on to a differential equation on P;^, with respect to the 
covering (/? : Pi|. — ^ P^ determined by the rational function ip{x). 

Pull-back transformations between hypergeometric and Heun equations give 
identities between the classical Gauss hypergeometric and Heun functions. For 
example, we have 



Hn 



9 
7/9 



1/3,1 
7/9; 2/3 



ip{x) , 



(5.7) 



where (p{x) is as in Example 



2.4 



8^2 



and 0{x) = (l-a;)~i/36 (l - x 
The transformation of singularities and local exponents for Fuchsian equations 
is explained in 22, Lemma 2.1]. The prefactor 6{x) shifts the local exponents 
at some points, but does not change the exponent difference anywhere. The 
rational function ip{x) multiplies the local exponents and their differences 
by the branching order at each point. If Q is a singularity of the starting 
Fuchsian equation in d/dz, a point P above Q will be non-singular for the 
pulled-back equation only if the branching order at P is n and the exponent 
difference at Q is equal to 1/n (and Q is not a logarithmic point when n = 1). 
For example, the 2F1 function in (5.7) solves a hypergeometric equation in 
z,d/dz with the exponent differences 1/9, 1/2, 1/3 at 2; = 0, l,oo, respectively, 
while the exponent differences for the pulled-back Heun equation are the 
branching fractions 2/9, 1/3, 1/9, 2/3 at a; = 0, 1, t, 00, respectively. The roots of 
8a;'^ — 72a;^ — 27a;-|-27 became non-singular after the proper choice of d{x). By the 



1/12 



way, the rational function (p{x) of Example 2.4 is identified in our classification 
by the label B7 in Table 2.3.9 of f|3l]and in Appendix SJB] 

Recently, parametric transformations between Heun and hypergeometric 
equations without Liouvillian solutions]^ were classified in [52], [13]. They 
apply to hypergeometric equations where at least one exponent difference is not 
restricted to a value 1/n with n G N; hence a parameter. In total, there are 61 
parametric transformations up to the well known symmetries of hypergeometric 
and Heun equations [23]. But the number of Galois orbits of utilized Belyi 
coverings (up to Mobius transformations) is 48. These Belyi functions are listed 
in [HI Table 4]. They satisfy condition (i) but not (ii) of Definition 1.3 because 



the parameter(s) could be specialized to satisfy the hyperbolic condition. The 
parametric transformations are labeled P1-P61 in [23], following similar criteria 



^^Liouvillian solutions 124! of second order linear differential equations are the "elementary" 
solutions: power, algebraic, exponential, trigonometric functions, their integrals (in particular, 
logarithmic and inverse trigonometric functions). They can be written in the form y = 
exp(J r) 2Fi(</')) where r, ip are rational functions, 2F1 is Gauss' hypergeometric function with 
a reducible, dihedral or finite monodromy. There are algorithms to find Liouvillian solutions 
in this form [18], hence a table of pull-back transformations is not needed. 

The hyperbolic restriction l/A; + l/^+l/m < 1 gives a finite list of (fc, £, m)-minus-4-regular 
Belyi functions, while 1/fc + l/f+l/m > 1 would lead to infinitely many Belyi functions 
because of applications to Liouvillian functions. 
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as in Appendix here. The Belyi functions of this article complete the list 
of hypergeometric-to-Heun transformations when no Liouvillian solutions are 
involved. 

Remark 5.1. Non-existence of Belyi functions with some branching patterns 
can be proved by non-existence of implied transformations of Fuchsian equations 
§5/. For example, there is no (2, 3, 10)-mmMs-4 Belyi function with the 
branching pattern 9 [2] =6 [3] = [10] + 2 + 2 + 2 + 2, because it would pull-back 
a hypergeometric equation with the exponent differences 1/2,1/3,1/2 to a non- 
existent Fuchsian equation with a single singularity where the exponent difference 
is not 1 (but 5). This example illustrates that there are no (k, i,m)-minus-l 
functions, unless 1 S {k,£,m}. In the exceptions, the implied hypergeometric 
equation must have a logarithmic singularity with the exponent difference 1. In 
particular, the polynomial {x'^ + 1)''' is a {k,l,dk)-minus-l Belyi function, and 
{x^ — ixY^ /{x"^ — is a {2k, l,3k)-minus-l Belyi function. 

An interesting observation is that the pull-back covering (^(a;) can be 
recovered from local solutions of the related hypergeometric and Heun equations, 
if only an oracle would tell us one constant. Particularly, suppose that 
the point x = of Hcun's equation lies above the singularity z = of 
hypergeometric equation. Let yi , 2/2 denote the hypergeometric local solutions 



in (5.3), respectively, and let Yi,l2 denote the Heun local solutions in (5.4) 
(5.5), respectively. We have the formula Yi{x) — 9{x)yi{ip{x)) like (5.7) 



and a similar formula [32 Lemma 3.1] relating ^2,^2 but normalized by a 
constant K that depends on the first power series term of (p{x). The quotient 
iIji{x) = Y2/Y1 does not depend on the prefactor 9{x), and can be identified 
with the respective quotient ipoiz) = j/2/yi up to the constant multiple K. We 
have ^pi{x) = x^'" (1 + 0{x)) and Vo(^) = z^-^ {I + 0{z)). The identification 
ii^o{z) = Kipi^x) gives z — tpQ^^Ktpi^x)). Therefore, the Belyi covering ip{x) 
the composition of the inverse function ip^^ with "01(2;) multiplied by the proper 
constant K. For instance, the Belyi function of Example |2.4| can be computed 
by inverting the function 



^1/9^^/5/36,17/36 



and composing with 
{Kxy^^Rn 



V 10/9 



9 

187/81 



5/9,11/9 
11/9; 2/3 



/^p^ / 1/36, 13/36 



Hn 



8/9 



7/9 



1/3,1 
7/9; 2/3 



where K — —64/3. The ratio of two independent solutions of the same 
differential equation of order 2 is called a Schwarz map of the differential 
equation. We consider the Schwarz mapa^ again in Appendix §D| 



^^In the general context of Fuchsian equations related by a pull-back transformation, the 
pull-back covering can be similarly recovered by a proper identification (up to a constant 
multiple) of Schwarz maps as well. In fact, our implemented algorithms often assume a pull- 
back of a hypergeometric equation to a Fuchsian equation with 4 singularities (rather than 
canonically normalized Heun's equation), so to avoid extensions of the moduli field. This is 
done when two or more branching fractions are equal and represent points in the same fiber, 
as demonstrated by the polynomial W in Example |6.2| Instead of the constants t, q, K in 
96.11 the constants j, q, K were generally used. 
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This observation is significant in a few ways. Firstly, a data base of our Belyi 
functions could be given by the data of Heun equations to which they apply (the 
exponent differences, the parameters q, t), the hyperbolic type (fc, £, m), and the 
constant K. The Belyi coverings would be then recovered by reconstructing a 
rational function from a power series. If d is the degree of a Belyi covering, 
2d + 8 power series terms would suffice (and exclude most of false rational 
reconstructions). Secondly, given a branching pattern (and thus the exponent 
differences of presumably related Heun and hypergeometric equations), the Belyi 
coverings (p{x) could be computed by assuming undetermined constants t,q,K 
and finding algebraic restrictions between them for reconstruction of f{x) from 
the power series of 'ipQ^{Kipi{x)). This approach does not appear practical, 
but |6.2| presents a deterministic algorithm that uses an implied Heun-to- 
hypergeometric transformation in a similarly general way, and eliminates all 
undetermined variables except 3 before calling Grobner basis routines. And 
thirdly, our probabilistic algorithm §6.1| searches through all possible t, q, K 
in finite fields, reconstructs possible minus-4-hyperbolic Belyi functions over 
considered finite fields, and uses a version of Henzel lifting to produce Belyi 
functions in Q(x). 



6 The algorithms 

The list of minus-4-hyperbolic Belyi functions was originally generated by a 
probabilistic algorithm by a thorough examination of Heun functions and their 
Schwarz maps over some finite fields, an d lift ing, identifying the obtained Belyi 



functions in Q(a;). This is explained in [6.1 The complete list was generated 
by considering at most 7 finite fields Fp = Z/(p) for p < 960, though eventually 
we kept the algorithm running for total 100 primes. In principle, this does not 
ensure completeness of the list however. 



The deterministic algorithm in [6.2 takes a branching pattern as an input. 



and produces the Belyi coverings with that branching pattern. By using the im- 
plied Heun-to-hypergeometric transformations, smaller degree algebraic systems 
for undetermined coefficients are obtained than with straightforward methods, 
and with far less parasitic solutions |llj . The deterministic algorithm produced 
the same Belyi maps (up to Mobius transformations) as the probabilistic one. 
Completeness of our results is proved assuming correct implementation of the 
deterministic algorithm. 

As a practical matter of confidence, the completeness of results is foremost 
verified by the same output of the two independent algorithms. In addition, 
we did a combinatorial search to find all minus-4-hyperbolic dessins d'enfant, 
up to degree 36. This gives a verification of a large part (« 95%) of relevant 
branching patterns, covering sa 91% of obtained dessins. We also compared the 
list of Belyi functions with the r-field in M with Felikson's list [7 of Coxeter 
decompositions in the hyperbolic plane; see Appendix S|Dj This provides enough 
confidence in completeness of our results. 
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6.1 A probabilistic modular method 

The used probabilistic algorithm is based on the expectation that a Belyi 
function with any realization field will be properly defined over a p-adic field Qp 
for some prime p among a sequence of considered subsequent or random primes. 
Concretely, suppose that a Belyi function ip{x) pulls-back a hypergeometric 
equation to Heun's equation with specific parameters t, q, and that respective 
Schwarz maps of both equations are identified by a constant K as described after 



Remark 5.1 1ft, q, K are elements of a number field Q(a), then ip{x) G Q^{a){x). 
By Chebotarev's theorem p4], the minimal polynomial for a has a root in Fp 
for a positive density of primes p. The density is at least 1/ D, where D is the 
degree of the number field Q(q;). For all but finitely many of those primes, we 
will have a € Qp and t,q,K ^ Zp (the p-adic integers). The Belyi function (p{x) 
can be found as follows: 

( j) Consider all possible values t,q,C G Fp of the reduction of t, q, K modulo 
an (eventually) suitable prime p; 

{a) Reconstruct (p{x) as a rational function in Fp[a;] by identifying the Schwarz 
maps as described after Remark |5.1| We need the first 2d + 8 terms in 
the Schwarz maps "00(2;) S'^d i)i{x) to be in Fp, so p has to be sufficiently 
large. For example, if a local exponent difference is 1/3, then we need 
p > 3(2(i + 8) to ensure that local hypergeometric solutions have the first 
2d + 8 coefficients in Zp. For degree 60 coverings, the starting prime was 
907 > 7(2 • 60 + 8). 

{Hi) Use Hensel lifting to obtain an expression of ^p{x) in Qp; 

(w) Use LLL techniques to compute minimal polynomials of its coefficients, 
thus reconstructing (f{x) as an element of Q(a)(a;). 

Our strategy is as follows. Foe each branching pattern of Tables 2.3.7-3.4.4, 
we run through a sequence of primes p and the possible reduced values t,q,C G 
Fp[a;]. For each of the 0{p^) pairs of t, q we have to compute series expansions 
for the solutions of Hq and Hi. This is done rapidly using linear recurrences 
for coefficients of these solutions; Maple has the command gfun[difFeqtorec] for 
getting the recurrences. We expect ip to be in Fp[[a::]] for suitable primes p. If 
ipQ o (p matches Kipi in Qp[[a;]], then this poses certain necessary conditiontj^ 
on the p-adic valuations of the coefficients of ipo and ipi . We compute the series 
solutions of Hq and Hi to enough precision so that we can test these necessary 
conditions. This way, many pairs t, q can be discarded, and we typically end up 
with O(p^) pairs. Thus, the rational reconstruction step (ii) "only" needs to be 
called for 0{p^) combinations oit,q,K. 

If we find a ip mod p, we store it in a file. Another program will Hensel 
lift it, apply LLL reconstruction to Q{a){x), and compare with the already 
computed data base. Each Belyi map (p has a density S^p of suitable primes. 
The expectation number of times that the same (p will be found is then 100 • S^p. 
Unless the density is tiny, the likelihood that (p will be found is very high. 



-'-•'if (f = Xx"^ + - ■ ■ e 'Lp[[x\\ is substituted into i/iq = x'^'^ {I + aix + a2X^ + . . .), -with S Qp, 
and if the first ^ Zp is a„, and if ipi = A'^ia;'*i'"(l + b-^x + b-^x^ + ■ ■ .) is the result of the 
substitution, then the first bi ^ I 
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The smallest (5^ encountered was 1/6, for the H10-H14 covering^^ with the 
realization field Q(C7)- Most of the table was found after just two primes. The 
first 10 primes took about a week on Maple, running on 12 CPU's. Among the 
100 primes, each Belyi function was found at least 16 times. 

The modular method is quite slow, because 0{p'^) combinations of i,q,C 
have to be inspected for each p. But its advantage is low requirement of 
computer memory. This means that the computation can continue for weeks on 
end, without a risk that the computation will halt due to memory problems, and 
without human intervention (this is important, because if human intervention is 
needed in any of the steps, then, in a table with hundreds of cases, a gap would 
become likely). 



6.2 A deterministic algorithm 

A (A;, ^, m)-minus-4 Belyi function is determined by a polynomial identity 

P'' U + Q'^ V ^ R"" W, (6.1) 

where P,Q,R are monic polynomials in C[x] whose roots are the regular 
branchings, and U, V, W are polynomials whose roots are exceptional points 
with correct multiplicities. The Belyi function is then expressed as 

^M^y_!L l-in(x) = —. (6.2) 

The polynomials P,Q,R should not have multiple roots; W may be monic. 
The 4 exceptional points may be immediately set to a; = 0, l,i,oo, but that 
may introduce an extension of the moduli field if there are exceptional points 
in the same fiber with the same branching order. If there is a point above 
(fi € {0, 1, oo} that has a unique (possibly regular) branching order for that 
fiber, it is convenient to set that point to a: = oo without extending the moduli 



field. The degrees of the polynomials in (6.1 ) are determined by the branching 
pattern and the assignment of a; = oo. 

The most straightforward computational method is to assume undetermined 



coefhcients of the polynomials in (6.1), and solve the resulting system of 



algebraic equations between the coefficients. This is not practical for Belyi 
functions of degree > 12, mainly because of numerous parasitic [IV solutions 



where some polynomials in (6.1) have common roots. 

A more restrictive set of equations for undetermined coefficients can be 
obtained by differentiating ip{x), as in [531 §4.1]. In particular, the roots of 
f'{x) include the branching points above (f = 1 with the multiplicities reduced 
by 1. A factorized shape of logarithmic derivatives must be the following: 

^=h,^, 4^=h,^. (6.3) 
(fix) QRF ' ip{x) - 1 PRF ^ ' 

Here hi, /12 are constants, and F is the product of irreducible factors UV W, 
each to the power 1. On the other hand, 

^=e^ + ^~rn^~^ (64) 



-"^^The estimate 5^ > 1/deg Q(q) is sharp when Q{a) D Q is a Galois extension. This is the 
case for Q(C7). Higher degree encountered number fields (such as for J28) had significantly 
higher 5^ > 1/6. 
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and similarly for the logarithmic derivative of (p{x) — 1. Both logarithmic 
derivatives have two expressions. This gives a generally strongei[^ set of 
algebraic equations. If fc = 2, the polynomial P can be even eliminated 
symbolically. This setting does not use the location (p = 1 oi the first fiber. 
For convenience, we assume the polynomials U, V, W to be monic therefore. 



Then (p{x) in (6.2) is determined up to a constant multiple, which can be found 
by at the latest stage by evaluating (6.11 at a root of Q or V. If a; = oo is 
in the (p = oo fiber additionally, the constants /ii,/i2 in (6.3) are equal to the 
branching order at a; = oo, which is then the residue value of (6.4) at a; = oo. 



To get an even more restrictive system of algebraic equations, we utilize 
the fact that our Belyi functions transform hypergeometric equations to Heun 
equations. The method bluntly uses the following lemma. 



Lemma 6.1. Let 'p{x) he a Belyi map as in (6.2). Hypergeometric equation 



( 5.1 ) with 



.-In 



k e m 2 



1 1 

k^ £ 



1 

m 



1 



is transformed to the following differential equation under the pull-back trans- 
formation z I— > ip{x), y{z) H- > 



d''Y{x) 



dx^ 
'A 



B 



F' U' V W 

Y ^ kU ^ JV ^ mW 
hih2 P^-'^Q^-'^UV 



Y{x) 
dx 



R^F^ 
mR'W 
RW 



mR' 



2^,2 




\ mR" 
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kU 


IV 


W J 



Proof. A lengthy symbolic computation, using (6.2) and (6.4) 



Y{x) = 0. 



□ 



The transformed equation is to be identified with the target Heun equation, 
or (if the roots of U, V, W are not normalized to a; = 0, 1, t, oo) with a Fuchsian 
equation with 4 singularities at the roots of UVW . The accessory parameter q 
is a new undetermined variable. The terms to dY{x)/dx are always identical, 
but comparison of the terms to Y[x) gives new algebraic equations between 
the undetermined variables. If fc = 2, ^ = 3, not only P but also Q can be 



eliminated symbolically. The two expressions in ( 6.4 ) and Lemma 6.1 then give a 
non-linear differential equation for i?, with q and the coefficients of (presumably 
monic) /7, V, W as parametric variables. After substitution of general polynomial 
expressions for R and C/, V, W, we collect to the powers of x and get a system 
of algebraic equation for undetermined coefficients. 

Example 6.2. Consider computation of degree 54 Belyi functions with the 
branching fractions 1/7,1/7,1/7,2/7. We assign the branching fraction 2/7 



^^Particularly, the new system has less parasitic solutions in general. Let G = 
gcd{P''U,Q''V,R'^W) and H = PQRF/G. Then only those parasitic solutions remain for 
which H is a. polynomial. Besides, a root of G is then a root of H if and only if it divides 
one of P^U,Q^V, R"^W in a higher order than others, and then it must be a simple root of 
H. This follows from invariance of the rational functions in | |6.3| in the hyperbolic and the 
reduced settings. 
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to X = oo, so that U = V ^ 1, and the polynomials P,Q, R,W are monic, 
without multiple roots, of degree 27, 18, 7, 3 respectively. If we would assume 
W — x{x — l){x — t), the Heun equation would have a — 9/14, b — 13/14 and 
c = d = 6/7. To avoid increase of the moduli field, we rather assume W = 
x^ + W2X + W3. Here the x^ term is zero-ed by a translation x ^ x + f3, so that 
only scaling Mobius transformations x ^ ax are left to act. The transformed 
Fuchsian equation must have the following term to Y{x): ab{x—q)/W. Symbolic 
elimination ofP,Q on Maple gives the following differential expression: 

7R"" 7R'" (W R'\ {lR"f R"W' ( IZW 35i?' 



IbR m \W R J 26i?2 RW V 7W 13R 

SR" 1033 \ /3,_„ . 16W''^\ nW 



115g - —X + —- -(163a: - 247g) 



7RW V 13 y 7R^W V2' ' 131^ J 2W^ 

3i?' //183 241 \ W 67\ 18 / 13 \ /46 27 



9 - — + - + -— 2a; - -g -X - -<z = 



2RW VV 7 13 J W 2lJ \ 5 7 Vl3 7 

Here the values W" — 6x, W'" — 6, W"" = are simplified. Substituting the 
explicit W , R = x^+ria;^ + . . .^r^x+r^i and clearing the denominator, we obtain 
a polynomial expression of degree 15 in x. The leading term gives q — — 5ri/52. 
The next term gives nothing nei^^ But the next 5 equations allow subsequent 
elimination of ^3, r4, rs, rg, ry in terms of ri, r2, u)2, W3. The 4 remaining 
variables are weighted-homogeneous, with weights 1,2,2,3. Elimination of 
W2,W3 using the other 10 equations is done with the Grobner basis routine of 
Maple 15 in about 35s. The algebraic system has 4 Galois orbits of solutions, 3 
of them parasiti^^ In the proper solution, we can take 

W^x'- 4899a; - 370078, R = x' + 28x'' + ^90632^5 ^5 + . . . . 

The expression for ip{x) is long. We looked for an optimizing Mobius 
transformation. The bit size of (p{x) is reduced by the factor ~ 2.26 after the 
Mobius substitution x i-> (241a; — 212)/(a; + 4). Then 

ip{x) 



86A{x - 4)(3a;2 + l)(a; + 4)2i?7 ' 

where i? = 3a;^ - 7a;^ - 14^^ - 98a;'* + 147x^ - 7x'^ + 56a; + 16 and 

Q = 47a;i* - 2028a;*^ + 55G2a;i^ + 54540x*^ - 263535a;*^ - 32592a;" + 2249268a;^^ 
-3436872a;" + UU5x^° - 1425900^*^ - 8774370a;** - 1715652a;^ - 10594017a;^ 
+2223144a;^ - 5284080a;'* + 1638144a;^ - 1306368x^ + 239616a: - 135168. 



^^We observed the following property of the polynomial system generated by just j6.4| l. Let 
El denote the sequence of polynomial equations given by the subsequent leading x-terms from 
the identification of ip' {x)/(p{x), and E2 denote the same sequence from the other logarithmic 
derivative. If x = 00 is in the fiber = 00 and has branching order e, then the first e equations 
from El and the first e — 1 equations from E2 are usually independent, but the eth equation 
from E2 is always a consequence of those. 

^'^The parasitic solutions are: the degree 18 covering mentioned in footnote 3; a degree 10 
covering with the branching pattern 5[2] = 3[2] + 4[1] = 7 + 2 + 1; and the non-cyclic cubic 
Belyi covering. In all cases, the simplification of the numerator and the denominator of ip{x) 
is by a linear polynomial to the maximal power (36, 44 or 51). 
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Even if {k, i) ^ (2, 3), most of the coefficients of P, Q, R can be eliminated 
subsequently from the two expressions of logarithmic derivatives and the Y(x) 
term of the transformed Fuchsian equation. If a; = oo is assigned as a point above 
If = oo, the highest order terms usually allow elimination of all but 3 variables 
(or all except 4 weighted homogeneous variables, if the scaling transformations 



X ^ ax are allowed to act). Lemma 6.1 fails to give new algebraic relations 
when A = 0. 

The order of k,£,m (or the 3 fibers) is not essential, of course. In our 
algorithm realization, we sought to assign a point a; = oo to an exceptional 
or regular point (preferably with a unique branching order in its fiber, and of 
maximal possible branching order e), and then assign that fiber as f = oo. With 



this we take advantage of the explicit constants hi = h2 = e in (6.4) and sooner 
eliminations. The hardest Grobner basis computation is with the last 3 (or 4 
if weighted homogeneous) variables. The strategy of using first a total degree, 
then elimination of 2 variables appears to be fastest for complicated examples. 
On Maple 15, the degree 60 Belyi maps were computed in 110s, the Galois orbit 
J28 in 274s, and the orbit pair Hll, J26 in 830s. 



Lemma 6.1 was already used with k — £ — 2 in [21] basically, to derive 
linear differential expressions for the polynomial components of Belyi maps 
for Kleinian transformations of hypergeometric equations with a finite dihedral 
monodromy group. 

6.3 Additional algorithms 

The main work for this paper was to compute all Belyi functions (up to Mobius 
equivalence, and up to Gal(Q/Q)-conjugacy) for the Heun table. But we 
computed a lot of additional data as well, and that required implementing 
numerous additional algorithms. The following is a summary, sorted roughly 
w.r.t. the amount of work involved. 

(1) Given a 4-hyperbolic Belyi map, compute its branching type, its i-value, 
j-invariant, the canonical realization field, and moduli field. 

(2) Given (51,50,500) and (5i,5o,5oo), with 5^,5., e Sd, 5o5iffoo = 1 and 
<5o,5i> transitive (and likewise for the gi). Decide if they represent the 
same dessin or not (decide if 3heSd ^^^di^ ~ 9i * ^ {0, 1, co}). 

(3) Gompute the obstruction conic and a conic- model (if it exists) . 

(4) Decompositions. Given (5i,5o,5oo) representing a dessin, compute the 
subfield lattice of Q{x)/Q{ip), where each subfield is represented by its 
dessin, by which is meant the following: If Q{f ) C L C Q{x) is a subfield, 
then L = Q(5) for some 5 by Liiroth's theorem, and (p — h{g) for some 
Belyi map h. By dessin of L we mean the dessin of h. 

(5) Given a Belyi map in / g K{x) and an embedding K C compute the 
dessin d'enfant of the image of / under this embedding. 

(6) Given (p G K{x), if there exists a Mobius-equivalent if of substantially 
smaller bitsize, then find such ip. 
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Item (1) is done directly from the definition. Most of the time for item (1) 
was spent on optimizing the polynomials representing the number fields (this 
can be done with the polredabs command from GP/PARI, but we had to work 
around some efficiency problems, avoiding unnecessary integer factorization). 

For item (2), suppose that h{l) — b for some yet to be determined 
be{l,...,d}. Then -3^(1)) =ffo^ffr---5r''(&) for aU . . . , n,.. 
That determines h because <go,gi> acts transitively. So we can find h (if it 
exists) by checking d cases, b = 1,2, ... ,d. 

Item (3) was explained in Section [4] For item (4) we computed the 
subgroups H of the monodromy group G :~< go, gi > that contain {g G 
G\g{l) — 1}. Given such H, writing down the action of gi, go, g^o on the 
cosets of H produces the dessin of the subfield corresponding to H. We then 
identified the component Belyi maps (such as h of item (4)) by using the full 
table of known hypergeometric and hypergeometric-to-Hcun transformations; 
This way we obtained all decompositions of all entries of our table; see 
[171 Decomposition_or_GaloisGroup] for used notation and detailed composition 
lattices. 

Although algorithms for computing the monodromy exist |19j item (5) 
was still a considerable amount of work because we had to develop our own 
implementation, specifically optimized for rational functions of high degree that 
ramify over only 3 points. 

Item (6) consists of a collection of algorithms, we will sketch a few. Suppose 
ip G K{x) and S = (p~^{{0, l,oo}). If some a G S has a minimal polynomial 
/ G K[x] of degree 3, we can apply polredabs to find an optimized polynomial 
g, then compute a Mobius transformation over K that will send a to a root of 
g, and check if it makes ip smaller. A similar trick will work if we have three 
points of degree 1, or one of degree 1 and one of degree 2. 

Given a factor / = a„a;" + • • • apx^ we can also try scaling, i.e. multiply x 
by primes that appear in a„ , oq to see if that makes ip smaller (implementing 
this for number fields is work) . Scaling is generally most effective after applying 
other steps first. For example, if S has just one point in ii'lJIcx)}, then first 
shift that to oo, then select a factor / = a„a:" + • • • aox'^, clear its second highest 
term with x t-^ x — a„_i/(na„), and then apply scaling. 



A Appendix: Sorting criteria 

In |3.1[ the minus-4-hyperbolic Belyi functions were grouped into 10 classes 
A-J. We order the Belyi functions inside those classes by the following criteria: 

(a) the first criterium is the j-invariant; 

(b) the second criterium is the branching fractionj^ 

(c) the last criterium is the degree of the covering. 

The sort of j-invariants lexicographically adheres to the following criteria: 
(al) the j-field; 

^*The first two criteria establish that our hst is basically sorted by Heun equations. To 
identify the Heun equations, invariants describing accessory parameters should be added |23l 
§D]. 
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(a2) the i-field; 

(aS) the leading coefficient of the minimal polynomial in Z[a::] for the j- 
invariant. 

The order of j-fields and t-fields is settled by the following criteria: 
(fl ) the field degree; 

(f2) if the field is a quadratic extension of Q then: 

(}la) real quadratic fields have precedence over Q(-\/a) with a < 0; 
(flh) the fields Q(-\/a) with the same sign of a are ordered by the increasing 
order of \a\. 

(fS) if the field is of higher degree, then the criterium is the field discriminant. 
The integers in (aS) and (fS) are ordered as follows: 

(11) the product of the primes dividing the integer; 

(12) the absolute value. 

The numbers in (il), (i2), (fib) and (c) are ordered in the increasing order. The 
tuples of branching fractions are ordered as follows 

(bl ) in each tuple, the four branching fractions are ordered in the increasing 
order of their denominators, then secondarily the numerators. 

(b2) the tuples are compared lexicographically, from their first elements, and 
the elements are matched first by their denominators then numerators. 

These criteria break all ties in our list of Bclyi functions. Due to (il), the fields 
or i-values that ramify or degenerate modulo the same set of primes are placed 
next to each other. The leading coefficient in (aS) gives information about the 
primes where the covering is ramified. In particular, for j G Q the leading 
coefficient is the denominator of j. 

B Appendix: The A-J tables 

The following pages display tables of Galois orbits of minus-4-hyperbolic Belyi 
functions, grouped as specified in §3.1| and ordered by the criteria in j fA) All 
tables contain the following columns: 

• Id: the label from Al to J18; 

• Branching frac: the branching fractions of a minus-4-hyperbolic function; 

• d: the degree of a Belyi function; 

• [Mm]: the values of k,£,m written compactly. For k — 2,£ — 3,m > 10, 
only the value of m e [10, 14] is given. 

• Monodromy/comp. or Mndr/cmp.: The monodromy group G — ... 
is given for indecomposable coverings, and compositions are indicated 
otherwise. The composition notation is explained in ^Cj 
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T J 

Id 


Brandling irac. 


a 


[kim\ 


Monodromy/comp. 






Al 


1/2, 1/2, V5, V5 


6 


[255] 


G = A 






— 


-1 


A2 




12 


[245] 


2 [445] o 6, Al o 2 




— 


-1 


A3 


1/2, 1/2, 2/7, 2/7 


18 


[237] 


2 [277] o 9 






— 


-7 


A4 


1/2, 1/4, 1/4, 2/5 


12 


[245] 


2 [445] o 6 






-1 


-1 


A5 


1/2, 1/4, 1/8, 1/8 


24 


[238] 


2 [488] 2 [288] o 2 [248] o 3 


— 





A6 


1/3, 1/3, 1/3, 1/9 


16 


[239] 


4 [339] o 4 






— 





A7 


1/3, 1/3, 2/3, 2/7 


16 


[237] 


2 [337] 8 






-3 


-3 


A8 


1/3, 1/3, 1/4, 1/4 


10 


[334] 








— 


-2 


A9 




20 


[238] 


2 [388] o 10, A8o2 


— 


-2 


AlO 


1/3, 1/3, 1/7, 3/7 


32 


[237] 


4 [337] o 8 






— 


-3 


All 


1/3, 1/3, 2/7, 4/7 


20 


[237] 


2 [377"] o 10 






— 





A12 


2/3, 2/3, 1/7, 1/7 


16 


[237] 


2 [337] o 8 






— 


-3 


A13 


2/3, 1/4, 1/8, 1/8 


20 


[238] 


2 [388] o 10 






-2 


-2 


A14 


2/3, 1/7, 1/7, 4/7 


20 


[237] 


2 [377"] o 10 












A15 


2/3,2/7,2/7,2/7 


20 


[237] 


2 [377"] o 10 












A16 


1/4, 1/4, 1/4, 1/4 


8 


[248] 


2h o 2 [444] o 


2 







-1 


A17 




12 


[334] 


2h o 2 [444] o 


3 







3 


A18 




24 


[238] 


see diagram ( 


C.l 


1 





— 


A19 


1/5, 1/5, 1/5, 1/5 


12 


[255] 


2h o Al 










A20 




24 


[245] 


2h o A2{A19, [445]} 






A21 


1/7,1/7,2/7,2/7 


48 


[237] 


2 [777] o 3 [337] o 


% 






A22 


1/7,1/7,2/7,4/7 


36 


[237] 


2 [777"] o 2 [277] c 


9 


-7 


-7 


A23 


2/7,2/7,2/7,2/7 


36 


[237] 


4 {A3, [277"] > 


}[277] o9 




-7 


A24 


1/9,1/9,2/9,2/9 


24 


[239] 


2 [999] o 3 [339] o 4 




-3 




33 



Id 


j-invariant 


Branching frac. 


d 


[Mm] 


Monodromy/comp. 


Bl 




1/2, 1/4, 1/8, 3/8 


18 


[238] 


3 [288] 2 [248] 3 


B2 




1/2,1/7,2/7,3/7 


27 


[237] 


3 [277] 9 


B3 




1/4, 1/4, 1/4, 1/4 


6 


[344] 


2^ 3 


B4 






12 


[246] 


4if{D25, P22{B3}, P39}o3 


B5 


2^73^/3* 


1/3, 1/3, 1/6, 1/6 


12 


[246] 


6{P16, [366]}[266] o2 


B6 




1/3, 2/3, 1/7, 2/7 


24 


[237] 


3 [337] 8 


B7 




1/3,2/3,1/9,2/9 


12 


[239] 


3 [339] 4 


B8 




1/4, 1/4, 1/5,4/5 


10 


[245] 


G ~ Sio 


B9 




1/4, 1/4, 1/8, 1/8 


10 


[248] 


G = Aq:C2 


BIO 






30 


[238] 


3 [388] 10, B9 3 


Bll 




1/6, 1/6,1/12,1/12 


18 


[12] 


6{P17, [3 12 12]} [2 6 12] o3 


B12 




1/7, 1/7,3/7,3/7 


36 


[237] 


4 [277] 9 


B13 




1/7,2/7,2/7,4/7 


30 


[237] 


3 [377"] 10 


B14 




1/8,1/8,3/8,3/8 


24 


[238] 


8{D9, [288]}[248] o3 


B15 


6481^/3^5^ 


1/2, 1/4, 1/4, 3/8 


15 


[238] 


5 [248] 3 


B16 




1/3, 1/7, 2/7, 4/7 


28 


[237] 


G = A2S 


B17 




1/3, 1/7, 3/7, 5/7 


16 


[237] 


G = Ai6 


B18 


73l273/223'^5^ 


1/3, 1/3, 1/3, 1/3 


5 


[335] 


G = Ar, 


B19 






10 


[10] 


2ho5, B18o2 


B20 


7^2287^/2*^325*^ 


1/3, 2/3, 2/7, 3/7 


12 


[237] 


G = 5*12 


B21 




1/4, 1/4, 1/4, 1/4 


6 


[344] 


G = S5 


B22 






12 


[246] 


2ho C24, B21 2 


B23 




1/5,1/5,2/5,3/5 


12 


[245] 


6 [255] 2 


B24 




1/5,1/5,1/10,3/10 


18 


[10] 


6 [2 5 lOJ 3 


B25 


2 7 ol 271 /o 11 

4993 /2 3 7 


1/3,2/3,1/7,4/7 


12 


[237] 


G = S12 


B26 


1/4, 3/4, 1/5, 2/5 


8 


[245] 


G — 5*8 


"DOT 




Vo> /O) /d> /o 


Q 



[ZDOJ 


lt = FbUd, Z ) : C/2 










[9/1 fil 


9 n R97 n 9 


B29 


9403761 no3/c6ii 4 


2/3 1/4 i/s 3/8 


14 


[238] 




B30 


2*181^2521^3654^^34 


1/3,1/4,1/8,5/8 


16 


[238] 


G = Aie 


B31 


TS'' 193^409^ 2^ 3^ 5'' 7*^ 


1/7,2/7,3/7,5/7 


18 


[237] 


G = ^is 


B32 


4920172^3*^52114 


2/3,1/7, 1/7,4/7 


20 


[237] 


G = S20 


B33 


2* 106791301^31452^8 -1^6 


1/2,1/7,3/7,4/7 


15 


[237] 


G = Si5 


B34 


829^30469^365678 ]^94 


1/3,1/7,2/7,5/7 


22 


[237] 


G = S22 



34 



Bl 



B2 






B3 



B4 



B5 



B6 



(XXX) 

B9 BIO 



B7 B8 




C7 




Bll 



B13 



oooo 



B14 



B17 



B21 




B15 



B16 






B18 , 



B19, 



B20 



B22 





B23 




B24 




B25 



B29 




B26 B27 

o CX>- 




B28 t-^; <3 



B30 




B31 




B32 




B33 




B34 



35 



Id 




j-invariant 


Branching frac. 


d 


r 7 /] 1 

[kim\ 


Monodroniy / comp. 


CI 


2 


2-^3-^ll-^ 


1/2, V2, 1/4, 1/4 


12 


[238] 


2ij 2 [248] 3 


C2 




5^ 11 ■''31 '72'"' 7*^' 


1/14,1/14,1/14,1/14 


18 


[14] 


2ho9 


C3 




7336013/2-3^7^ 


1/2, 1/2, 1/6, 1/6 


8 


[246] 


G =PSL(3,2):C2 


C4 






1/7, 1/7,3/7,6/7 


18 


[237] 


G = Sis 


C5 


3 


221933/3 


1/4,1/4,3/8,3/8 


18 


[238] 


2ho3 [248] 3 


C6 




2'^533/33 


1/3, 1/3, 2/7, 2/7 


32 


[237] 


4 [337] 8 


C7 






1/3, 1/3, 2/9, 2/9 


16 


[239] 


4 [339] 4 


C8 




33531573/2^11^ 


1/2,1/8,3/8,3/8 


15 


[238] 


G = 5'i5 


C9 




13'^541'^/3'^11 


2/3, 1/9, 2/9, 2/9 


14 


[239] 


G = *S'i4 


CIO 




109391333 /2**3^13'' 


1/4,1/4,1/5,3/5 


14 


[245] 


G — Si4 


Cll 






1/4, 1/4, 1/6, 1/6 


14 


[246] 


G =PSL(2, 13): (72 


C12 




333731926373/11^17'* 


2/3, 1/7, 1/7, 3/7 


26 


[237] 


G = 826 


C13 




2^" 5^1301^43889^317]^ ]^i0j^34 


1/3,1/11,2/11,2/11 


16 


[11] 


G = Aie 


C14 


5 


2" 


1/2, 1/2, 1/4, 1/4 


10 


[245] 




C15 




2^173 


1/2,1/5,1/5,1/10 


15 


[10] 


5 [2 5 10] 3 


C16 






1/2,1/5,2/5,2/5 


10 


[245] 


5 [255] 2 


C17 






1/4, 1/4, 1/4, 1/4 


20 


[245] 


4h{F13x,C14}o5 


C18 




1036813/3^5 


1/7, 2/7,3/7,3/7 


30 


[237] 


G = S30 


C19 


6 


2^97l3/3^ 


1/2,1/4, 1/4,3/5 


8 


[245] 


G = As 


C20 




2^1934673/3^5^ 


1/3,2/3,1/4, 1/4 


6 


[334] 


G = Se 


C21 








12 


[238] 


C20o 2 


C22 




11312593/2-335'' 


1/3, 1/3, 1/8, 1/8 


26 


[238] 


G =PSL(2,25):C2 


C23 






1/9,1/9,2/9,5/9 


18 


[239] 


G = Sis 


C24 




11319793/233-512 


1/2, 1/2, 1/4, 1/4 


6 


[246] 


G = S5 


C25 






1/10, 1/10,3/10,3/10 


18 


[10] 


2ho9 


C26 


7 


3353173 


2/3,2/7,2/7,3/7 


14 


[237] 


G = Si4 


C27 




2'*37327l3/3'55'' 


2/3, 1/7, 1/7, 5/7 


14 


[237] 


G = *S'i4 


C28 




22^^31073/3127 


1/3,1/3,3/8,3/8 


14 


[238] 


2ho7 


C29 




2'^5'^16073/3i'^7^ 


1/3,2/3,2/3,1/7 


8 


[237] 


G^As 


C30 


10 


79493/2^31° 


1/2, 1/2, 1/4, 1/4 


10 


[245] 


G^Aq 


C31 






1/6, 1/6, 1/6, 1/6 


10 


[256] 


2h 5 


C32 




II3133233/2-3125 


2/3,2/3,1/8,1/8 


10 


[238] 


2ho5 


C33 


13 


1122973/2^32013 


1/2,1/3,3/7,3/7 


13 


[237] 


G = Ai3 


C34 


21 


331273/5S 


2/3,1/7,3/7,3/7 


14 


[237] 


G — Si4 


C35 




3313673/2^52 


1/2,2/7,3/7,3/7 


15 


[237] 


G = S'15 


C36 




7573ll8273/2"*3'^17*' 
37 5653 // 3 5 7 


1/2,1/7,3/7,3/7 


21 


[237] 


G = A21 


C37 




1/3,1/3,1/5,1/5 


1-7 

7 


food 

[335] 


G = Ai 


C38 








14 


[10] 


2ho7, C37o2 


C39 


105 


243^86813/5^ 


1/2,2/7,2/7,3/7 


21 


[237] 


A21 


C40 




33273293/2^5-7^ 


2/3,2/3,1/7,2/7 


10 


[237] 


G = Sio 


C41 


273 


5^349-''8516873-y2-''"3''' 7" 11*513 


1/2,1/3,2/9,2/9 


13 


[239] 


G = Ai3 


C42 


385 


332891893/2185711 


1/2,2/3,2/7,2/7 


11 


[237] 


G = Sii 



36 



CI C2 



C3 C4 



C5 



CIO 



C18 



C22 



C31 





C7 C8 C9 

0>o- o^>- 00<] 



Cll 







1 ""'^ 





C13 



C16 C17 



C14 CIS C16 




C26 C27 C28 C29 



C32 C33 




C34 





C39 ^-^ C40 C41 C42 

0*0 K>0- 00- 




37 



Id 


V 


j-invariant 


Branching frac. 


d 


\}ilm\ 


Monodromy/comp. 


Dl 


-1 


223^1375^ 


1/4, 1/4,1/8,5/8 


18 


[238] 


6 [248] 3 


D2 




-2*1093/5*^ 


1/3, 2/3, 1/9, 1/9 


14 


[239] 


G — tS'i4 


D3 






1/5,1/5,1/5,4/5 


12 


[245] 


6 [255] 2 


D4 






1/5,2/5,1/10,1/10 


18 


[10] 


6 [2 5 10] 3 


D5 






1/7, 1/7, 2/7,5/7 
1/10, 1/10, 1/10, 1/10 


30 


[237] 


G — 5*30 


D6 






24 


[10] 


4//{G5x,P43}o6 


D7 


-2 


-25193/3*^ 


1/3, 1/3, 1/3, 1/3 


8 


[334] 


2ff o4 


D8 








16 


[238] 


4^/{G13x,P41{D7}}o4 


D9 




2-47V3^ 


1/2,1/2,1/8,3/8 
1/6, 1/6, 1/6, 1/6 


12 


[238] 


4 [248] 3 


DIO 






8 


[266] 


2ff o4 


Dll 








16 


[246] 


4jf{Gllx,P23{D10}}o4 


D12 




-2*^2393/31° 


1/2, 1/3, 1/4, 1/4 

1/3,1/3,1/3,2/3 
1/3,1/6,1/12,1/12 


8 


[246] 


4 [344] 2 


D13 






8 


[238] 


4 [334] 2 


D14 






16 


[12] 


4 [3 4 12] 4 


D15 




-48264l3/2-32llio 


1/11,1/11,2/11,3/11 


18 


[11] 


G — S\% 


D16 




-2549773/2^31^194 


1/3,1/3,1/8,3/8 


20 


[238] 


G = A20 


D17 




7607^1753^2 ''3205"'ll"' 


1/3,1/7,1/7,4/7 
1/2,1/6,1/6,1/6 


34 


237 


G — 


D18 


-3 





12 


[246] 


3 [366] 2 [2661 o 2 

L J L J 


D19 






1/3, 1/3, 1/3, 1/3 


12 


[239] 


3 [339] o 4 


D20 






1/3,1/3,1/3,3/7 

1/4,1/12,1/12,1/12 
1/7, 1/7, 1/7,6/7 
2/7, 2/7, 2/7, 3/7 
1/8,1/8,1/8,3/8 


24 


[237] 


3 [337] o 8 


D21 






18 


[12] 


3 [312 12] o 2 [2612] o 3 


D22 






30 


[237] 


3 [377"] o 10 


D23 






30 


[237] 


3 [377"] o 10 


D24 






30 


[238] 


3 [388] o 10 


D25 




2II/3 


1/2,1/2,1/4, 1/4 


6 


[246] 


2//o3 


D26 




-2173373/134 
-233713/3 • 78 


1/3,1/3,1/7,6/7 
1/2, 1/4, 2/5, 2/5 
1/7, 1/7,1/7,2/7 

1/7, 3/7, 3/7, 4/7 


14 


[237] 


G = Si4 


D27 




9 


[245] 


G = A9 


D28 






54 


[237] 


G = 554 


D29 






18 


[237] 


G = Sis 


D30 


-5 


-528I3/3I65 


1/3,1/3,2/3,3/7 
1/3, 1/3, 1/4, 1/4 


10 


[237] 


G = Sio 


D31 






5 


[344] 


G = S5 


D32 








10 


[246] 


2ff o5, D31o2 


D33 




27914233/3<'577^ 

-113888113/26345-712 


1/3, 2/7, 2/7, 5/7 


16 


[237] 


G = Ai6 


D34 




1/4, 1/4, 1/7, 2/7 


10 


[247] 


G = Sio 


D35 




-11^23830621091^2*3^° 5^7* 43* 


1/3,1/7,1/7, 1/7 


52 


[237] 


G = A52 


D36 


-6 


-2^6359^2999^3751674 


1/5,2/5,1/6,1/6 
1/2,1/5,1/5,2/5 


8 


[256] 


G = Ss 


D37 


-7 


-53I6373/2I87 


7 


[255] 


G = S7 


D38 








14 


[245] 


D37o2 


D39 






1/3,1/3,1/3,2/9 

1/7, 1/7,2/7,4/7 


14 


[239] 


G = Si4 


D40 




-53373l673/283'ill'i 


36 


[237] 


G = Ase 


D41 




-2653144II3/36731110 


1/3,1/11,1/11,3/11 


16 


[11] 


G = A16 


D42 


-14 


-2^199287631^/3265673 


1/3,2/5,1/10,1/10 


16 


[10] 


G = 


D43 


-15 


-2693/21035 


1/2,1/5,1/5,3/5 


10 


[245] 


5 [255] 2 


D44 






1/2,1/10,1/10,3/10 
1/4, 1/4,1/4,1/4 


15 


[10] 


5 [25 10] 3 


D45 






20 


[245] 


2h C30 


D46 




-113593/2123-53 


1/2,1/7,1/7,4/7 


27 


[237] 


G = 1527 


D47 




-333350893/2145^234 


2/3,1/7,2/7,2/7 


26 


[237] 


G = 526 


D48 


-35 


1685104151'y26332577. ^^34 


1/3,1/3,1/10,3/10 


14 


[10] 


G = Su 


D49 


-39 


-17^29^5197^2303352133 


1/2,1/9.1/9,4/9 


15 


[239] 


G = Si5 


D50 






1/2,1/10,1/10,3/10 


15 


[10] 


G = Si5 



38 



Id 


discQ(i) 


j-invariant 


Branching frac. 


d 


[klm] 


Monodromy/ comp. 


El 


-2^3^ 


-2633233 


1/2, 1/3, 1/3, 1/3 


6 


[334] 




E2 








12 


[238] 


El 2 


E3 




-33173/2 


2/3, 1/9, 1/9, 1/9 


18 


[239] 


G = Sis 


E4 


-2^3* 


-33533833/2^ 


1/3, 1/3, 1/3, 1/4 


9 


[334] 


G=((C3XC3):Qs):C3 


E5 








18 


[238] 


E4o 2 


E6 






1/3,1/3,1/3,4/7 


18 


[237] 


G = Sis 


E7 


-2-^5^ 


-5^^2413/23 


1/2,1/10,1/10,1/10 


18 


[10] 


G = Sis 


E8 






1/5,1/5,1/5,2/5 


10 


[255] 


G = Sio 


E9 








20 


[245] 


E8 2 


ElO 




5-2113/215 


1/2, 1/2, 1/2, 1/4 


5 


[245] 


G = S5 


Ell 






1/4, 1/4, 1/4, 1/4 


20 


[245] 


4ij{H17xx} oElO 


E12 




-2^33613/5^ 


1/3, 1/3, 1/3, 1/5 


6 


[335] 


G = A6 


E13 








12 


[10] 


E12o2 


E14 






1/3, 1/3, 1/3, 5/7 


12 


[237] 


G = Ai2 


E15 


-2*^11'^ 


-2*112133/36 


2/3,1/11,1/11,1/11 


14 


[11] 


G = Si4 


E16 


-2^13^ 


-33413833/2-137 
-2^5*113/35 


1/13, 1/13, 1/13, 2/13 


18 


[13] 


G = Sis 


E17 


-2^335^ 


1/5,1/5,1/5,3/5 


8 


[255] 


G = As 


E18 








16 


[245] 


E17o 2 


E19 


-2*337* 


-2i"'7i 173233/313 

03-72 

—2-3 7 


1/3, 1/3, 1/3, 2/3 


14 


[237] 


G = Si4 


E2U 


o9o3'74 

—2 3 7 


V3, 1/7, 1/7, 1/7 


1 A 

10 


[237] 


(_T — 


E21 


-263773 


-3-2233/28 


1/2,1/2,1/2,2/7 


9 


[237] 


G = Sg 


E22 






2/7 2/7 2/7 2/7 


36 


237 


4h{H37xx} oE21 


E23 


-2^5373 


2^3337393/5" 75 


1/5, 1/5, 1/5, 1/7 


8 


[257] 


G = As 


E24 


-265374 


337-20993/21*57 


1/2, 1/4, 1/4, 1/4 


7 


[247] 


G = St 


E25 


-265^^133 


-5-3410909y2^''3i''l3^ 


1/2, 1/4, 1/4, 1/4 


15 


[245] 


G = 5i5 


El 




E2 E3 ^ 


-T— ^ E4 




E5 






E6 



€>0 



E7 




E8 






E9 




ElO o 



Ell 



E16 



E12 



E13 



E14 



E15 




E21 



40 



T J 

Id 




aisc<Q[t) 


Brancnmg irac. 


7 

a 


r 7 /7 1 
[Ktml 


Monodromy/comp. 


Fl 


2 


as j 


1/3, 1/3, V3, V3 


8 


[334] 


G =PSL(3,2) 


F2 








16 


[238] 


2h F5, Fl 2 


F3 






1/3,2/7, 3/7, V7 


16 


[237] 


G = 


F4 




2" 


1/2, 1/2, 1/8, 1/8 


18 


[238] 


G =PSL(2,17) 


F5 




1 

2"3^ 


1/2, 1/2, 1/3, 1/3 


8 


[238] 


G =PSL(3,2):C2 


F6 




2^17 


1/8, 1/8, 1/8, 1/8 


36 


[238] 


2h F4 


F7 




-2*^327 


1/2, 1/3, 1/3, 1/4 


7 


[334] 


G =PSL(3,2) 


F8 








14 


[238] 


F7o 2 


F9 


3 


as j 


1/3,1/3,1/12,1/12 


14 


[12] 


G ==PSL(2,13):G2 


FIO 


5 


as j 


1/2,1/2,1/10,1/10 


12 


[10] 


G =PSL(2,11):C2 


Fll 






1/5, 1/5, 1/5, 1/5 


12 


[255] 


G =PSL(2,11) 


F12 








24 


[245] 


2//OF16, Fllo2 


F13 




-2^52 


1/2, 1/2, 1/4, 1/4 


10 


[245] 


2ho5 


F14 




-2*5^ 


2/3,1/5,1/10,1/10 


14 


[10] 


G — Si4 


F15 






1/4, 3/4, 1/5, 1/5 


12 


[245] 


G^Ai2 


F16 




2*^5^ 


1/2, 1/2, 1/5, 1/5 


12 


[245] 


G =PSL(2,11):G2 


F17 




2*5211 


1/10, 1/10, 1/10, 1/10 


24 


[10] 


2h FIO 


F18 


6 


-2*^3^5 


1/3,1/3,2/3,1/4 


5 


[334] 


G^S's 


F19 








10 


[238] 


F18o2 


F20 


7 


2*3^^72 


2/7, 2/7, 3/7, 3/7 


24 


[237] 


2h F22 


F21 




-263-72 


1/3,1/3,1/7,3/7 


32 


[237] 


G = A32 


F22 




-283-72 


1/2,1/2,2/7,3/7 


12 


[237] 


G ~ S12 


F23 


21 


as j 


2/3, 1/7, 2/7, 3/7 


20 


[237] 


G = S20 


F24 




-3372 


1/7, 1/7,3/7,5/7 


24 


[237] 


G = A2i 


F25 


22 


2^3-113 


1/2,1/3,1/3,3/8 


11 


[238] 


G = Sii 




Id ^ discQ(t) Branching frac. d \klm\ Monodromy/comp. 



Gl 
G2 
G3 
G4 
G5 
G6 
G7 
G8 
G9 
GIO 
Gil 
G12 
G13 
G14 
G15 
G16 
G17 
G18 
G19 
G20 
G21 
G22 
G23 
G24 
G25 
G26 
G27 
G28 
G29 
G30 
G31 
G32 
G33 
G34 
G35 
G36 
G37 
G38 
G39 
G40 
G41 
G42 
G43 
G44 
G45 
G46 
G47 
G48 
G49 
G50 
G51 
G52 



-2 



-3 



as J 

28 
2^5 



2^5 

2^13 

2^5-13 

as j 

2" 

2103 

2"3 

263-11 
2^3-19 

12:2='^3'' 

as j 

2^33 

337 

2^327 



337-13 

12: 3^7^13'' 
243-537 

as j 

23^2 
3472 
2273 

2572 

223273 
227211 



7237 

12:2=7i°43^ 
-11 12:3''11» 

-15 2^3^52 



223353 



72 

1/3 

2/3 

1/4 



1/2 
1/2 
1/3 
1/2 
1/3 
1/2 
1/3 
1/2 
1/2 
1/8 
1/2 
1/2 
1/8 
1/3 
1/4 
1/3 
1/3 
1/3 
1/3 
1/2 
1/2 
1/2 
1/3 
1/2 
1/2 
2/3 

2/3 

1/7 
1/3 
1/2 
1/3 
1/2 
1/2 
1/2 
1/7 
1/2 
1/3 
1/2 
1/2 
1/3 
1/2 
1/5 



/2 



253352 Y4 

12: 2^3i»53l93 1/2 



172 

1/5 
1/4 
1/4 



1/8 
1/2 
1/4 
1/8 

2/3 

1/2 

2/3 

1/2 
1/2 
1/8 
1/8 
1/2 
1/8 
1/7 

3/4 

1/3 

3/4 

1/3 
1/6 
1/7 
1/3 
1/12 
1/13 

3/4 

1/7 
1/7 
1/7 
1/7 
1/7 
1/2 
1/3 
1/4 
1/2 
1/2 
1/7 
1/7 
2/7 
1/7 
1/7 
1/3 
1/2 
1/5 

1/2 

1/4 
1/9 



7V78 15 [238r 
16 [10] 



1/10, 3/10 
1/8, 1/8 
1/8,1/8 

1/8, 1/8 
1/10, 1/10 
1/8, 1/8 
1/8,5/8 
1/8,3/8 
1/6, 1/6 
1/8, 1/8 

1/3.1/3 
1/8,3/8 

3/8,3/8 

1/8,3/8 
1/2,1/8 
1/8, 1/8 

2/7,6/7 

1/6,1/6 
1/9,4/9 
1/8, 1/8 
1/12, 1/12 
1/12, 1/12 
1/7,6/7 
1/3, 1/3 
1/12,1/12 
1/13,1/13 
1/5, 1/5 

2/7,4/7 
2/7,4/7 

1/7,2/7 
1/7,4/7 
2/7, 2/7 
1/2, 1/3 
1/3, 1/3 
1/4, 1/7 
1/7,4/7 
1/3, 1/3 

4/7,4/7 

2/7, 2/7 

2/7,4/7 

1/7,2/7 
1/7, 1/7 
1/3,1/11 
1/2,1/5 
1/5, 1/5 

1/5,1/5 

1/5,2/5 
V9, V9 



20 [238] 

10 [248] 

30 [238] 

27 [238] 
12 [10] 

28 [238] 
15 [238] 
12 
8 



[238] 
[246] 
18 [238] 



8 
12 
24 



14 
16 



9 
15 
16 

7 
21 
14 
32 
42 
40 

7 

28 
8 
12 
14 
24 
33 
22 
39 
45 
12 
6 

12 
24 

6 



[235 
[235 
[238] 



21 [238] 
9 [238] 
36 [238] 
16 [237] 
8 [246] 

14 [239] 
16 [238] 

[12] 
[12] 

15 [237] 
[239] 
[12] 
[13] 
245 
237 
237 
237 
237 
237 
237 
237 
247 
237 
237 
237 
237 
237 
237 
237 

[11] 
245" 
255 
245' 

[255^ 



12 [245" 
18 [245^ 
21 [239' 



5 [248[ o 3 

G = ^16 
G = 520 

2ho5 

2// o Gl: 5 :G4 o 3 

G = S'27 
2ff 06 

G = ^28 

G = S'15 
G = S12 
2ff o4 
G = Si8 
2^0 4 
G — S12 
2h o G14 

G = ^21 

G = {{{C3f:Qs):C3):C2 
4h{J1xx}oG17 

G = A16 
G = 5*8 
G = Si4 
G = A,e 
2ff o7 
G = A16 
G — S'15 

G=PSL(2,8):G3 
G = S'15 

G = A16 
G = S7 

G = ^21 

G = Si4 

G = S32 

G = S42 

G = A40 

G =PSL(3,2) 

4ff{I3x,G39}oG35 

G = {C-if :PSL(3,2) 

G = S12 

2h o G35 

2jj oG38 

G = A33 

G = S22 
G = S39 

G = A45 
G = M12 
G^Ae 
2h o G49 

4H{I11^ G50{G47}}oG46 
G = A 

2h o G46, G49 o 2 
G = Si8 
G = A2i 



42 



Id discQ(j) discQ(i) Branching frac. d [k£m\ Monodromy /comp 

xJi o3 „ ,, „■ I/O I/O 1 /fi 1 /A ori room r)OT /o i n\ 



3^ 



G =PSL(2,19) 

G = S12 
G = Sio 
G = S'15 

G={A6:C2):C2 
G = A,e 

G = 5*12 
2jj oH7 

G = ^22 

G =PSL(2,29) 
G =PSL(2,43) 
G =PSL(2,13) 
2h o H12 
2hoH10 
G — Ss 

2//OH19, H15o2 
2h o ElO 

G={{C2)^:As):C2 
G = Ss 

G = Aio 



HI 
H2 
H3 
H4 
H5 
H6 
H7 
H8 
HQ 
HIO 
Hll 
H12 
H13 
H14 
H15 
H16 
H17 
H18 
H19 
H20 
H21 
H22 
H23 
H24 
H25 
H26 
H27 
H28 
H29 
H30 
H31 
H32 
H33 
H34 
H35 
H36 
H37 
H38 
H39 
H40 
H41 
H42 
H43 
H44 
H45 
H46 
H47 
H48 
H49 
H50 
H51 
H52 
H53 



as J 
as j 

-223105 



2233 


as j 




-28365 


2-^ 3"* 


as j 


2^3* 


-2«3« 




2^3^11 




-2^:iHi 


2 


2675 




3375 




263374 




7^13 




7*29 


2^52 


as j 




2"54 




21254 




-21454 




-21455 




293355 




-28335^^7 


3^5 


-36537 


2211 


-2^^11313 




-2^3^7-113 


2^13 


as j 


2^3-5"^ 


-2^3^5'^ 


22345 


-24385311 


3^11 


-3^113 


527 


5572 




-265572 


223-72 


-24335-7'' 


223^7 


283672 




283773 


283.72 


-283374 


225-7 


as j 




285273 




-24335372 




-24537219 




24537319 


225-72 


18._230gl07l2 


2^5-72 


-2I33-5274 




283-5375 


7-172 


-3-5-72175 


227-13 


as j 




-243-72133 


3-7-11 


as j 


223-7211 


as j 



1/3 
1/3 
1/3 
1/2 
1/2 
1/3 
1/2 
1/9 
1/3 
1/2 
1/3 
1/2 
1/3 
1/7 
1/5 

1/2 
1/2 
1/2 
1/2 
1/3 

1/3 
1/2 
1/3 
1/2 
1/2 
1/3 
1/2 
1/3 
1/2 
1/2 
1/3 
1/2 
1/2 
1/3 
1/2 
1/3 
1/2 
1/2 
1/4 
1/3 
1/2 
1/2 
1/2 
1/7 
1/2 
2/3 
1/2 
1/2 
1/3 
1/3 
1/2, 



1/3, 

2/3 

1/3, 
1/3, 
1/4 
1/9 
1/2, 
1/9, 
V9 
1/2 
1/3 
1/2, 
1/3, 
1/7, 
1/5, 

1/2, 
1/4, 
1/2, 
1/2, 
1/3 

1/3, 
1/5, 
1/9 
1/11 
1/3, 
1/3, 
1/3, 
2/3, 
2/3, 
1/3. 
1/3 
1/2, 
1/3, 
1/7, 
1/2, 
1/3, 
1/2, 
1/7, 
1/4, 
1/3, 
1/3, 
1/7, 
1/2, 
1/7, 
1/2, 

2/3 

1/3, 
1/3, 
1/3 
1/7, 

2/3, 



/9, 1/9 
1/9,2/9 
2/3, 1/9 
1/9, 2/9 
1/4, 1/6 

2/9,4/9 

1/9,2/9 
2/9,2/9 
1/9,2/9 
1/7, 1/7 
1/7, 1/7 
1/3, 1/3 
1/3, 1/3 
1/7, 1/7 
2/5,2/5 

1/4, 1/4 
1/4, 1/5 
1/5,2/5 
1/3, 1/4 
1/4, 1/4 

1/8,5/8 
1/5, 1/6 
1/9,5/9 
1/11,2/11 

1/3,4/7 
1/11,2/11 

1/4,3/8 
1/10, 1/10 
1/9, 1/9 
1/3,2/9 

3/7,3/7 

1/3,3/7 
1/3,2/3 
1/7,6/7 
2/7,2/7 
2/3,2/7 
1/4,2/5 

2/7,5/7 
2/5,2/5 

1/7,5/7 
2/7. 2/7 
1/7,5/7 
1/2, 1/7 
1/7, 1/7 
2/3, 1/7 
1/7, 1/7 
1/3,3/7 

2/7,4/7 
2/7,3/7 

3/7,4/7 
V7,3/7 



20 
12 
10 
15 
10 
16 
12 
24 
22 
30 
44 
14 
28 
60 
8 

16 
10 
16 
8 

10 
10 
20 
14 
7 

16 
15 
11 
14 
13 
12 
11 
11 
20 
10 
7 

22 
18 
16 

7 
15 
14 
20 
25 
21 
15 
60 

8 

16 
17 
13 
26 
22 
11 



239 
239^ 
239 
239^ 
246' 
239' 
239 
239' 
239' 

23r 

237 
237 
237 
23r 

255' 
245" 
245' 
245' 
245' 
238' 
334' 
238 
238 
256' 
239 

[11] 
[237] 

[11] 
[238] 
[10] 
239 
239' 
237' 
237 
237 
237' 
237' 
237' 
245" 
237 
245 
237' 
237' 
237 
237 
237 
237 
237 
237 
237' 
237' 
237' 
237' 



G = 


Si4 


G = 


S7 


G = 


Aie 


G = 


S15 


G = 


Sii 


G = 




G = 


Al3 


G = 


S12 


G = 


Sii 


G = 


Sii 


2ho 


H34 


G = 


^10 


G = 


^7 


G = 


^'22 


2ho 


E21 


G = 


•^le 


G = 


^7 


G = 




2ho 


H39 


G = 


A20 


G = 


A25 


G = 


A21 


G = 


Ai5 


4jj{J19x 


G = 


Ss 


2ho 


H47 


G = 


AiT 


G = 


Al3 


G = 


S26 


G = 


S22 


G = 


Sii 



44 




45 



Id 




discQ(j) 


discQ(i) 


Branching frac. 


d 


[kim] 


Mndr./cmp. 


11 


4 


2'^3'' 


as j 


V2, 


1/3, 1/4, V6 


9 


[246] 


G = 


A9 


12 




22^3 


as j 


V3, 


2/3,1/7,2/7 


24 


[237] 


G = 


S24 


13 






2123276 


V2, 


1/2, 1/3, 1/3 


14 


[237] 


2ho 


G35 


14 




2873 


21376 


V2, 


1/2,1/7,2/7 


24 


[237] 


G = 


5*24 


15 








V7, 


1/7,2/7,2/7 


48 


[237] 


2ho 


14 


16 






24.23756^22 


2/3, 


1/7, 1/7, 1/7 


38 


[237] 


G = 


'5'38 


17 




2672 


2123.5.75 


V2, 


1/4, 1/7, 1/7 


9 


[247] 


G = 


A9 


18 




3'^7 


3*^527213 


V3, 


1/7, 1/7,5/7 


28 


[237] 


G = 


A28 


19 






24.3307l7;^g3 


V2, 


1/3, 1/3, 1/3 


21 


[237] 


G = 


A21 


110 




-2^11 


as j 


V2, 


2/3, 1/4, 1/8 


11 


[238] 


G = 


Sii 


111 






2i^3<^5^ 


V2, 


1/2, 1/5, 1/5 


12 


[245] 


2ho 


G46 


112 






as j 


V2, 


1/2, 1/4, 1/6 


7 


[246] 


G = 


S7 


113 






2123773 


V4, 


1/4, 1/6, 1/6 


14 


[246] 


2ho 


112 


114 




223373 


24:2^*3^^7^*13'* 


V3, 


1/3, 1/3, 1/7 


36 


[237] 


G = 


^36 


115 




-2^3^72 


as j 


V2, 


1/3,2/3,2/7 


9 


[237] 


G = 


5*9 


116 




2^3-11 


as j 


V3, 


1/4,1/8,3/8 


22 


[238] 


G = 


5*22 


117 




223'*13 


as j 


V2, 


1/3, 1/3, 1/9 


13 


[239] 


G = 


^13 


118 




2^3313 


as j 


V2, 


1/2, 1/3, 1/8 


13 


[238] 


G = 


^13 


119 






21436^33 


V3, 


1/3, 1/8, 1/8 


26 


[238] 


2ho 


118 


120 




7^19^ 


as j 


V2, 


1/3,2/7,3/7 


19 


[237] 


G = 


Sw 


121 




2^3^5-7-13 


as j 


V2, 


1/3,1/7,5/7 


13 


[237] 


G = 




122 


5 


2II73 


-2^1 76 


V2, 


1/2,1/3,2/7 


16 


[237] 


G = 


5*16 


123 






2263376 ;^3 


V3, 


1/3,2/7,2/7 


32 


[237] 


2ho 


122 


124 






as j 


V2, 


1/4,1/5,3/5 


9 


[245] 


G = 


^9 


125 




2^3^ 7 


as j 


V2, 


1/7,2/7,3/7 


27 


[237] 


G = 


5*27 


126 




223572 


as j 


V3, 


2/3,1/7,3/7 


18 


[237] 


G = 


5*18 


127 




335-73 


-2i°3'^527'^ 


V2, 


1/2, 1/7, 3/7 


18 


[237] 


G = 


^18 


128 






385376 ]^7 


V7, 


1/7,3/7,3/7 


36 


[237] 


2ho 


127 


129 




2^53112 


as j 


V2, 


1/2, 1/4, 1/5 


11 


[245] 


G = 


^11 


130 






2153.561^5 


V4, 


1/4, 1/5, 1/5 


22 


[245] 


2ho 


129 


131 




265^132 


as j 


V2, 


1/4,1/5,2/5 


13 


[245] 


G = 


Sis, 


132 




24327311 


-2*3''5^7''ll^31 


V2, 


1/7, 1/7,3/7 


33 


[237] 


G = 


^33 


133 




223.73^72 


as j 


V2, 


2/3,1/7,2/7 


17 


[237] 


G = 


517 




46 



Id 




discQ(j) 


discQ(t) 


Branching frac. 


d 


[Mm] 


Mndr./cmp. 


Jl 


6 




23^3^ 


1/2, 1/2, 1/8, 1/8 


18 


[238] 


2ho 


G17 


J2 




-2"3*'5^ 


as j 


1/2,1/3,2/3,1/8 


9 


[238] 


G = 


^9 


J3 




-3*5^7 


233I65573 


1/3, 1/3, 1/9, 1/9 


20 


[239] 


2ho 


J4 


J4 






2103175572 


1/2, 1/2, 1/3, 1/9 


10 


[239] 


G = 


^10 


J5 




-3-''5''^13 


as j 


1/2,1/3,1/5,1/10 


13 


[10] 


G = 


^13 


J6 




-2''335-7^ 


2^23752710 


1/3, 2/3, 1/7, 1/7 


30 


[237] 


G = 


Ssa 


J7 




-2''7''23 


21278233 


1/3, 1/7, 1/7, 2/7 


46 


[237] 


G = 


5*46 


J8 




-5^7^11 


3551174^^^3 


1/2,1/3,1/3,1/10 


11 


[10] 


G = 


^11 


J9 




-2''3'*7-ll'' 


283*7^11^13 


1/2,1/3,1/11,1/11 


13 


[11] 


G = 


^13 


JIO 




-223372172 


as j 


1/3, 1/7, 2/7, 3/7 


34 


[237] 


G = 


•534 


Jll 




2^3''7'^192 


as j 


1/2,1/3,1/7,4/7 


19 


[237] 


G = 


Sl9 


J12 




-2^7211213 


2^367^11^133 


1/3,1/3,1/7,4/7 


26 


[237] 


G = 


S26 


J13 


7 


-227'"^ 192 


-2^3*^7^°19^ 


1/3, 1/3, 1/7, 2/7 


38 


[237] 


G = 


S38 


J14 




-227'*233 


-2'*3^7^23^ 


1/2, 1/3, 1/3, 2/7 


23 


[237] 


G = 


S23 


J15 




-223'' 7^11 


-2^3i'55-7i°ll3 


1/3, 1/3, 2/3, 1/7 


22 


[237] 


G = 


S22 


J16 




-3^527*17 


-3''5^7*173l9 


1/3, 1/7, 1/7, 3/7 


40 


[237] 


G = 


^40 


J17 




-2''5-7"*17-23 


2854711^73233 


1/2, 2/3, 1/7, 1/7 


23 


[237] 


G = 


'S'23 


J18 


8 


2"5'^17 


2285"13-173 


1/2, 1/4, 1/5, 1/5 


17 


[245] 


G = 


Air 


J19 


9 


2125676 


236510714 


1/2, 1/2, 1/7, 1/7 


30 


[237] 


2ho 


H45 


J20 




22^132172 


-2^338134 ;^75 


1/2, 1/3, 1/3, 1/8 


17 


[238] 


G = 


Air 


J21 


10 


_2233452]^g2 


as j 


1/2, 1/3, 1/4, 1/8 


19 


[238] 


G = 


Sl9 


J22 




-2^3^5^7^1l2 


as j 


1/2,1/3,2/3,1/7 


15 


[237] 


G = 


Sl5 


J23 


11 


— 2''3^5^7° 


as j 


1/2, 1/3, 1/7, 3/7 


25 


[237] 


G = 


A25 


J24 


13 


2103579314 


as j 


1/2,1/3,1/7,2/7 


31 


[237] 


G = 


S3I 


J25 




223^7^11^19 


-2283127I811II192 


1/2, 1/2, 1/3, 1/7 


22 


[237] 


G = 


A22 


J26 






263197201111132193 


1/3, 1/3, 1/7, 1/7 


44 


[237] 


2ho 


J25 


J27 


14 


_28527l0ig.294 


2^*3"5*7^°13- 19^29^ 


1/2, 1/3, 1/3, 1/7 


29 


[237] 


G = 


^29 


J28 


15 


-2^^3^7^°11^31-37^ 




1/2, 1/3, 1/7, 1/7 


37 


[237] 


G = 


-437 
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Other occasional columns: 

• j-invariant: given if it is in Q \ {1728}, in a factorized form; 

• dj: the degree of the j-field (in tables I, J); 

• disc Q(j), disc Q{t): the field discriminants. If the extension Q{t) D Q{j) 
is of degree 6, the degree of the the t-Reld is indicated in the disc Q{t) 
column in a small underlined font. 

• ^ : indicates the quadratic extension of either the i-field (in Tables C, D) 
or of the j-field (in Tables F, G); 

• m-^: the quadratic extension for the moduli field (only in table A); 

• : the quadratic extension for the r-field (only in table A) . 

The tables are supplemented by pictures of respective minus-4-hyperbolic 
dessins d'enfant. We put « next to a few F, H-labels to indicate that their 
dessins are obtained by bi-coloring the vertices (and thereby making the edges 
thin). The J-pictures marked by the symbol # represent 4 dessins each, 
obtainable by reflecting (with respect to a horizontal axis) their left and right 
parts independently. The dessins of the following Galois orbits are displayed 
elsewhere: 

• H14, H46: in Figure [ij 

• H45, J19: in Figure [2) 

• B12, C6, C30, D45, F4, F6, HI, H10-H13: in Figure [3| 

• Fl, F7, Fll: in Figure|4l 



C Appendix: Composite Belyi functions 

Composition of a Belyi function (p{x) into smaller degree rational functions can 
be decided from the function field lattice between C(x) and C{lp), as described 



in [121 Section 1.7.2] and ^6.3 here. The subfield lattices are listed in our online 
table [TTl Decomposition_or_GaloisGroup]. 

On the other hand, composite minus-4-hyperbolic Belyi functions in- 
duce composite hypergeometric-to-Heun transformations. Thereby special 
cases of the parametric transformations P1-P61 of [23l §2.2] and the Heun- 
to-Heun transformations 2h, of [23l §4.3] often occur as composition 
parts. The quadratic transformation 2^ acts on the exponent differences as 
(1/2, 1/2, a, /?) ^ (a, a, /3, /3) and changes the j-invariant to a 2-isogenous j- 
invariant. The transformation Ah — 2h o 2h transforms (1/2,1/2,1/2,0;) ^ 
(a, a, a, a) and does not change the j-invariant. The composite transformations 
could be figured out by a careful consideration of possible compositions 
of hypergeometric-to-hypergeometric, indecomposable hypergeometric-to-Heun 
(parametric or some newly implied), and Heun-to-Heun transformations. That 



would constitute yet another checlj | of our list of Belyi functions 



^^In particular, any transformation to Heun's equation with 2 (or 3) exponent diflferences 
equal to 1/2 can be composed with 2h (or iu, respectively). Further, any Belyi function of 
the [fc^m]-type [344] or [266] gives rise to a type- [246] composition (with the degree doubled), 
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The most complicated decomposition lattice is for A18: 



= :P1 [288]:P2 



[4444] P12:[2244] [444] :P9 



CI: [2244] 



P4:[2224] 




248]^:A16 



[248]:P10 [238] 



334]:A17 



(C.l) 

In the square brackets, we see the [Mm] triples of intermediate hypergeo- 
metric equations, or similar indication of intermediate Heun equations. The 
transformation from [238] to the Heun equations is indicated before their 
square brackets. Similarly, the [k£m] triples are followed by the indication of 
transformations from them to the final [4444]. The x -power indicates two copies 
of that intermediate function field. The diagram includes PIO and P12, the most 
complicated parametric compositions [23,, §C]. The components [238]— [248] and 
[334] — [444] are cubic transformations, while the other lines represent quadratic 
ones (possibly 2h). 

In Tables of SjBj we indicate the components either by an A-J label from 
our list (if applicable), or by the degree otherwise. In the latter case, we 
give intermediate hypergeometric equations in the [k£m\ notation. Intermediate 
Heun equations are clear, hence no extras to 2^. Deeper branching is indicated 
by {}. The A-J, P labels inside them cither mean a transformation from a 
starting [fci?m] to an intermediate Heun equation (after Ah) or to the target Heun 
equation (otherwise). The xx -power indicates three copies of an intermediate 
function field. A label inside nested {} refers to a composition string avoiding 
the merging point of the outer {}. These hints should be enough to recover the 
composition lattices. 



D Appendix: Coxeter decompositions 



If a minus-4-hyperbolic Bclyi function in a canonical form (of Definition 2.2 1 
is defined over M, the Schwarz map^p"] of the related hypergeometric and 
Heun equations fit together nicely. Particularly, the quadrangle of Heun's 
equation is then tessellated into congruent (in the hyperbolic metric) triangles 



while all [334], [248]-type functions give type- [238] compositions, with the degree 2 or 3 times 
larger. In the same way, the [335], [255] Belyi functions give type [2 3 10], [245] (respectively) 
compositions. Quadratic transformation PI can be composed to CI and all compositions in 

Table A of fjB] as its j-invariant 1728 is 2-isogenous to itself and the j- value of CL 

■^''We alreaay considered Schwarz maps in the paragraph after Remark |5.l[ If a 
hypergeometric equation has real local exponent differences a,/?, 7 in the interval [0, 1], the 
image of the upper half plane C C under its Schwarz map is a curvilinear triangle with the 
angles 7ra, tt/?, 777. A nice illustration can be found in [S] pg. 38]. Analytic continuation 
of Schwarz maps follows the Schwarz reflection principle. Hodgkinson [9] first observed 
that pull-back transformations of hypergeometric equations induce tessellations of Schwarz 
triangles into smaller congruent Schwarz triangles. Similarly, if a Heun equation has real 
local exponent differences a,/3,7,5 in the interval [0,1], the image of its Schwarz map is a 
curvilinear quadrangle with the angles 7ra, 7r/3, 777, 7r(5. 
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of the hypergeometric equation. The degree formula in Lemma |3.1| ("m^ can 
be interpreted as the area ratio between the hyperbohc quadrangle and the 
triangles, if we multiply both the numerator and the denominator by tt. 

Subdivisions of hyperbolic quadrangles (or triangles) into congruent hyper- 
bolic triangles are called Coxeter decompositions in [7]. The list of Coxeter 
decompositions can be compared with our list of Belyi maps with the r-field C M, 
providing a mutual check of completeness. All Coxeter decompositions from 
our Belyi maps can be discerned in Figure [5] The similar pictures for Coxeter 
decompositions from parametric hypergeometric-to-Heun transformations are 
given in [22l Figure 2]. 

The Belyi functions of Tables D, E, G of ^|B]would certainly give no Coxeter 
decompositions, as their r-fields are not real. As shown in the last column of 
Table A, most of its Belyi functions have imaginary quadratic r-fields. Here is 
the listF^ of other extensions of the r- field over the i- field: 





B21, B22, Cll, C22, C24, F9 (and C6, F4, F6, HIO, H12); 




B9, BIO, D24, F16 (and C30); 




B6, B27, B28, C3, F2, F5 (and Fl, HI, Hll, H13); 




B2, C2 (and B12, H14); 




B18, B19, C31; 


-5+V5 


FIO, F12, F17 (and Fll); 



(and \/3 for D45). In the parentheses we list the obstructed Galois orbits of 
^ Evidently, the obstructed Belyi functions have no real r-fields. Among the 
non-mentioned Galois orbits of Belyi functions: 

• All other cases of the A, B, C-tables give a Coxeter decomposition each; 
10 + 23 + 34 in total. 

• Each entry of the F-table with discrim Q{t) < gives one Coxeter 
decomposition; 10 in total. 

• The entries F3, F23 give pairtj^of Coxeter decompositions as Q{t) = Q{j); 
F20 gives another pair with the i-field Q(\/7, ^3); and F25 gives none with 
the i-field Q(v/4^722^^). 

■^^The cases where the moduli field is an extension of the j'-field are: A4, A7, A13, A22 (as 
shown in Table A of ^ and B2, B6, D20, D24, F7. The coverings D20 and D24 have j = 0, 
but their moduli fields are 3) and (Q(\/— 2), respectively. The function F7 is considered 

in j ]4.5| It does give one Coxeter decomposition, into symmetric (tt/S, 7r/3, 7r/4) triangles. 

The list of cases with additional extensions for the r-field correlates well with the list of 
Belyi coverings with interesting monodromy groups (such as PSL in tables of ijB| and the list 
of multiple Galois orbits with the same branching pattern (as one can inspect empty entries 
in the first columns in tables of [13. 1|. 

^■^Dessins d'enfant and Coxeter decompositions are different geometric representations of 
a Belyi covering. The difference is twofold: the decompositions represent only a half of the 
Riemann sphere, and their vertices are the points not just above z S {0, 1} but above 2 = 00 
as well. To get a corresponding (real) dessin, two parallel copies of a Coxeter decomposition 
have to be glued along the edges to a topological sphere, and the vertices above z = 00, the 
incident edges and triangles have to be removed. A real conjugation in Q{t) D Q{j) acts on a 
Coxeter decomposition as a Mobius symmetry. 

Among the Galois orbits Al— J28 without an obstruction, the F-orbits and 110, 115, Jll 
have two real dessins, while other orbits have at most one real dessin. 



52 



• Each entry of the H-table with Q{t) — Q(j) or discrim Q(t) > gives a 
Coxeter decomposition; 11+11 in total. 

• Similarly, the odd degree I, J-orbits with Q{t) ~ Q{j) or discrim Q{t) > 
give Coxeter decompositions; 6+3 among 122-133 and 2 + 3 in the J-table. 

• The entries 110, 115, Jll with pairs of real dessins and Q{t) = Q{j) give 
pairs of Coxeter decompositions. 

In total, we must have 125 Coxeter decompositions, just as listed in f7\ Figures 
10 (5)~(11), 12, 13, 15-18]. There is a caveat, however. The decompositions 24 
and 36 in [7l Figure 18] coincide, while one triangulated quadrangle with the 
angles 7r/3, 27r/3, 7r/7, 37r/7 is missing. We identify the repeated decomposition 
as C4, and the missing one as 126. 

All Coxeter decompositions can be recognized in Figure [5j To assign Belyi 
functions (with the r-field in M) to Coxeter decompositions, we multiplying 
the branching fractions by n and look for quadrangles with the same angles. 
The pictures (a), (h) show the Coxeter decompositions 7, 6 in [71 Figure 15] 
into triangles with the angles 7r/2,7r/3 and 7r/12 or tt/II. They represent the 
Belyi functions Ell and C13, respectively. Picture (c) contains two hyperbolic 
quadrangles subdivided into twelve (7r/2, 7r/4, 7r/6)-triangles. They represent 
the Belyi functions B5 and B4, and coincide with the triangulations 3, 4 in [3 
Figure 12], respectively. Picture (d) contains the first five triangulations in [7l 
Figure 15], into (7r/2, 7r/3, 7r/10)-triangles. Here are the labels of Belyi maps 
and the quadrangles: 

C38: BCFK, F14: ACEG, C15: ACFH, C25: ACDG, B24: ACFG. 
Picture (e) contains the (7r/2, 7r/3, 7r/9)-triangulations in 7-, Figure 16]. Here is 
the respective sequence of Belyi maps and the quadrangles: 
H2: DFKM, B7: EGKM, C41: EHLM, C7: EHKM, C9: GFKM, 
H4: CFKN, H6: CFKP, A6: EFKM, C23: AGFK, H8: ABFK. 
There is initial ambiguity for assigning B7 and H2 because of the same 
branching fractions. But B7 is a composition 3 [339] o 4 as shown in the B- 
tablc, and its Coxeter decomposition splitf[^ into 3 triangles with the angles 
7r/3, 7r/3, 7r/9 (each formed by 4 smaller triangles). Further, picture (f) contains 
the (7r/2, 7r/3, 7r/8)-triangulations in T, Figure 17]. Here is the respective 
sequence of Belyi maps and the quadrangles: 



Bl: 


ABEK, 


A18: 


BEGK, 


B14: 


BERN, 


A5: 


BGKN, 


CI: 


ADST, 


110: 


FGKS, 


110: 


AGET, 


C32: 


CEKQ, 


F19: 


DFUQ, 


C21: 


ADFU, 


H29: 


ADFV, 


C28: 


ABRH, 


B29: 


AGEK, 


B15: 


ADSK, 


F8: 


ADFL, 


B30: 


ADFM, 


C8: 


ABET, 


H22: 


PDFM, 


C5: 


ADEK. 







Picture (g) contains the (7r/2, 7r/4, 7r/5)-triangulations in |71 Figure 13]. Here is 
the respective sequence of Belyi maps and the quadrangles: 



B26: 


FKPS, 


H39: 


VYQT, 


C19: 


VKQT, 


C16: 


AGOS, 


B23: 


AGQS, 


H16: 


AGQR, 


124: 


EHPS, 


H17: 


DHPZ, 


F13: 


VKPT, 


C14: 


BXZT, 


B8: 


BEST, 


H41: 


BFQT, 


129: 


OLNZ, 


F15: 


OLNS, 


CIO: 


VLNS, 


131: 


WUYQ, 


H18: 


WUKP, 


C17: 


BGPT, 


A20: 


AEMR, 


130: 


WULN. 



^•^ Coxeter decompositions do not always split according to (all) compositions of their Belyi 
functions, because smaller degree components do not necessarily have Coxeter decompositions. 
For example, consider A18 = A16 o 3, A19 = 2h o Al, B4 = 2^ o D25, B14 = 2^ o D9, 
J19 = 2h o H45, J26 = 2h o J25, etc. 
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There is the same ambiguity between C14 and H17 due to the same branching 
fractions. It is resolved by the composition C17 = 2//oC14. Picture (h) contains 
the (7r/2, 7r/3, 7r/7)-triangulations in [71 Figure 18]. Here is the respective 
sequence of Belyi maps and the quadrangles, with the repeated decomposition 
(36) replaced by the one for 126: 



C34 


CHRE, 


F3 


HMRYj, 


A15 


MTV A, 


All: 


HKRYs, 


H33 


HORY., 


A14 


A/STV, 


A21 


ELTW, 


C29 


REAA, 


115: 


MST.A, 


C42 


GMSY, 


F22 


CMsn, 


115 


UVAQ, 


B33 


CHRTL, 


C40: 


GARY., 


H53 


garh. 


B20 


mrj:a, 


BIT 


CHRX, 


B25 


CARX, 


H50: 


GMSX, 


C33 


HMSY, 


C26 


CMRT., 


C35 


CMRIi, 


C27 


GARY, 


C4: 


GHRY, 


H48 


AARX, 


F3 


CMRX, 


H40 


CMSY, 


133 


AMSX, 


H37: 


BMSY, 


B32 


HRYT, 


B31 


C'MRY, 


B34 


EHRX, 


F23 


AMRX, 


C39: 


BMRY, 


F20 


AMRY, 


126 


NRY^A, 


F23 


CNRY., 


C36 


GNRU, 


H52: 


GNRX, 


F24 


GNRY, 


F20 


CMTV, 


B13 


CMTW, 


123 


EKRX, 


Jll: 


DHRX, 


Jll 


HLSY, 


H44 


EHRX, 


125 


BNRY, 


J24 


FKRY, 


J19: 


FLSY, 


H36 


GHRY, 


C12 


ANRX, 


J26 


GMSY, 


B16 


GMRY, 


C18: 


ANRY, 


J17 


EQUW, 


J23 


GPQY, 


128 


ANRZ, 


F21 


GORY. 











The ambiguity between A14 and B32 is resolved by the reflection symmetry of 
A14 = 2o 10. The non-parametric decompositions (5)-(ll) of [71 Figure 10] and 
the decompositions (1), (2) of [71 Figure 12] represent the Galois orbits F18, 
B3, C20, C37, F7, H21, A17, H15, A19, respectively They can be obtained 
from our listed quadrangles of (respectively) F19, B4, C21, C38, F8, H22, A18, 
H16, A20 by pairing their triangles to larger triangles with the requisite angles 
(7r/3,7r/3,7r/4), (7r/3, 7r/4, 7r/4), (7r/3, 7r/3, ^5) or (7r/2, 7r/5, 7r/5). 



E Appendix: Arithmetic observations 

As observed in [531 §2.3], the t-parameters of Heun equations reducible to 
hypergeometric equations by a pul-back transformation are arithmetically 



interesting. The whole orbit (2.1) of values can be encoded by an arithmetic 
identity A + B = C with algebraic integers A, B, C (as "co-prime" as possible), 
as the set {A/C, B/C, C'/A, C/B, -A/B, -B/A}. Here are these identities for 
a few i-orbits in Q: 

B25 : 1 + 2-112 = 3^ B29 : 2^ 11^ = 5^ B30 : l-f 3^5^ = 2^13^ 
B31 : 1 + 2^3-52 = 7"^, B33 : 11^ + 2^7'' = 3^5, B34 : 7^^ + 3^5^ = 2^^192. 



The terms in these identities involve only small primes, usually in some 
power. Correspondingly, the ^-values factorize nicely in Q. These identities 
are interesting in the context of the ABC conjecture [21j and S-unit equations 
[M| . The "factorization" pattern holds for the t- values in algebraic extensions 
of Q as well, though arithmetic quality is then measured more technically [16] 
by the prime places and arithmetic height in p2(Q). The underlying reason is 
that the Belyi coverings (of pull-back transformations) tend to degenerate only 
modulo a few small primes P]. Hence the t-orbit (2.1) degenerates only modulo 
those bad primes. 

Amidst the encountered examples, we find the following well-known idcnti- 
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ties A + B = C in quadratic fields: 



C18:^ +2^3275 = 



2 

D37/D39 : f iif^j +7^='^"^"^ 

They are among top 12 known examples of remarkable ABC identities (TB] 
in algebraic number fields. Their ABC-quality is a; 1.697794, 1.707222, 
respectively, while Nitaj's table [TB] includes examples with the quality > 1.5. 
The Belyi function D42 gives a new example in Q(-\/— 14) with the quality 
log(3"53)/log(56 • 2 • 7 • 32 • 52) w 1.581910. However, the class number of 



-14) is equal to 4, hence an explicit arithmetic identity is less impressive, 
without 13th powers: 

(5 - 2^^) (11 + ^/^)' + {V^f = (5 + 2^^) (11 - ■ (E.l) 

Less symmetric quadratic identities arrise from the F, G-cases with Q{t) = Q{j). 
For example, G30 gives 

(i±^) " + (2 + v^f = 0. (E.2) 

The Belyi coverings ElO/Ellgive the following A + B ~ C example in a number 
field of degree 6. Let C denote a root of + 4z^ — 3z^ + 2. Then 



c + c _ 5c-^ + c y / i-c 3C" - 3c^ + c - C" 

2 4 / \ 2 4 



C + C' 5C^ + c'\^Vi + C 3C2 + 3C3 + C* + C'^ 



The numbers under the 23rd power have the norm 2, while the numbers in the 
(-6)th power are units. 
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